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FLIOHT-PATH  CHARACTERISTICS  FOR  DECELERATING 
FRCM  SUPERCIRCULAR  SPEED 
By  Roger  W.  Luidens 


SUMMARY 

Characteristics  of  the  following  six  flight  paths  for  decelerating 
from  a  supercircular  speed  are  developed  in  closed  form:  constant  angle 
of  attack,  constant  net  acceleration,  constant  altitude,  constant  free- 
stream  Reynolds  number,  and  "modulated  roll."  The  vehicles  were  required 
to  remain  in  or  near  the  atmosphere,  and  to  stay  within  the  aerodynamic 
capabilities  of  a  vehicle  with  a  maximum  lift-dreig  ratio  of  1.0  and 
within  a  maximum  net  acceleration  G  of  10  g's.  The  local  Reynolds 
number  for  all  the  flight  paths  for  a  vehicle  with  a  gross  weight  of 
10,000  pounds  and  a  60°  swept  wing  was  found  to  be  about  0.7X10®. 

With  the  assumption  of  a  laminar  boundary  layer,  the  heating  of  the 
vehicle  is  studied  as  a  function  of  type  of  flight  path,  initial  G  load, 
and  initial  velocity.  The  following  heating  parameters  were  considered: 
the  distribution  of  the  heating  rate  over  the  vehicle,  the  distribution 
of  the  heat  per  square  foot  over  the  vehicle,  and  the  total  heat  input 
to  the  vehicle.  The  constant  G  load  path  at  limiting  G  was  found  to 
give  the  lowest  total  heat  input  for  a  given  initial  velocity.  For  a 
vehicle  with  a  maximum  lift-drag  ratio  of  1.0  and  a  flight  path  with  a 
maximum  G  of  10  g's,  entry  velocities  of  twice  circular  appear  thermo¬ 
dynamically  feasible,  and  entries  at  velocities  of  2.0  times  circular 
are  aerodynamically  possible.  The  predominant  heating  (about  85  percent) 
occurs  at  the  leading  edge  of  the  vehicle.  The  total  ablated  weight  for 
a  10,000-pound-gross-weight  vehicle  decelerating  frcm  an  initial  velocity 
of  twice  circular  velocity  is  estimated  to  be  5  percent  of  gross  weight. 
Modifying  the  constant  G  load  flight  path  by  a  constant -angle-of -at tack 
segment  through  a  flight-  to  circular-velocity  ratio  of  1.0  gives  essen¬ 
tially  a  "point  landing"  capability  but  also  results  in  an  Increased 
total  heat  input  to  the  vehicle. 


INTRODUCTION 

Studies  of  interplanetary  missions  indicate  that  atmospheric  braking 
can  be  markedly  more  efficient  (in  terms  of  weight  requirements)  than 
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propulsive  braking  (ref.  l) .  To  apply  atmospheric  braking  at  the  entry 
velocities  associated  with  interplanetary  trips,  several  problems  must 
be  studied.  One  of  these  is  the  question  of  whether  the  achievable  aero¬ 
dynamic  entry  corridors  (e.g.,  refs.  2  and  3)  are  compatible  with  the 
corridor  capability  of  the  planet  approach  guidance  system.  A  second 
problem  is  protecting  the  vehicle  from  the  heat  generated  in  decelerating 
from  the  entry  velocity. 

In  berms  of  acceptable  G  loads,  entry  to  Earth  at  velocities 
greater  than  circular  appears  feasible.  The  advent  of  ablation  heat- 
protection  systems,  in  contrast  to  radiation-  and  heat-sink  systems,  also 
makes  possible  the  consideration  of  supercircular  entry  speeds  from  the 
thermodynamic  point  of  view.  To  minimize  the  weight  of  the  heat- 
protection  system  requires  the  use  of  "deceleration  flight  paths"  that 
have  inherently  low  total  heat  inputs.  The  deceleration  flight  path  is 
defined  as  that  segment  of  the  total  flight  path  In  the  atmosphere  which 
is  flown  for  the  purpose  of  decelerating  the  vehicle  from  an  initial 
velocity  and  altitude  to  a  desired  terminal  condition.  It  is  subsequent 
to  the  maneuver  which  determines  the  entry  corridor.  The  velocity  change 
during  the  entry  maneuver  may  be  small  (ref,  3),  so  that  to  a  first  ap¬ 
proximation  the  initial  velocity  for  the  deceleration  flight  is  equal  to 
the  entry  velocity. 

A  number  of  studies  of  the  heat  input  and  flight  mechanics  of  de¬ 
celeration  flight  paths  have  been  made  (refs.  4  to  /).  The  stagnation- 
point  heat  input  per  unit  area  nondimensionallzed  by  functions  of  nose 
radius  Rj,^,  lift-drag  ratio  L/D,  wing  loading  W/a,  and  drag  coefficient 

Cp  is  the  parameter  usually  presented.  Items  not  considered  in  these 

studies  which  are  dealt  with  in  the  present  analyses  are:  (l)  the  inter¬ 
relation  among  the  parameters  Rj^,  L/D,  W/A,  and  Cj)  required  by  the 

vehicle  geometry,  (2)  the  heat  input  to  the  windward  side  of  the  vehicle, 
(5)  the  integration  of  the  heat  input  per  unit  area  over  the  vehicle 
area,  and  (4)  the  effect  of  differentiating  between  the  vehicle  maximum 
and  operating  lift-drag  ratios.  The  parameter  heat  input  per  unit  ve¬ 
hicle  gross  weighs  is  used  as  a  measure  of  the  heat  input  to  the  vehicle. 
This  parameter  is  closely  related  to  the  ablated  weight  as  a  fraction  of 
the  vehicle  gross  weight. 

The  present  analysis  yields  closed-form  results.  The  flight  mechan¬ 
ical  approximations  made  are  less  restrictive  than  the  assumptions  used 
to  arrive  at  exist  ng  closed-form  solutions.  The  characteristics  of 
the  following  types  of  flight  paths  are  derived:  (l)  constant  angle  of 
attack,  (2)  const-  nt  G  load,  (3)  constant  heat  rate,  (4)  constant  alti¬ 
tude,  (5)  constant  free -stream  Reynolds  number,  and  (6)  "modulated  roll." 
The  following  constraints  are  imposed:  that  the  vehicle  (l)  remain  in 
or  near  ilie  atmr.sphere,  (2)  stay  within  Its  aerodynamic  capabilities, 
and  (w)  stay  wii.hin  the  G  load  capabilities  of  the  payload  and 
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structure.  A  Reynolds  number  constraint  is  also  discussed.  The  total 
heat  Input  to  the  vehicle  and  the  G  load  history  are  discussed  as  a 
function  of  the  type  of  flight  path,  initial  velocity,  and  the  initial 
G  load.  Finally,  a  method  of  controlling  landing  location  is  discussed. 


ANALYSIS 

The  general  objective  of  this  analysis  is  to  develop  closed-form 
solutions  for  the  characteristics  of  the  deceleration  flight  path. 


Equations  of  Motion 

The  equations  of  motion  Including  gravitation  and  aerodynamic  forces 
In  the  directions  tangent  and  normal  to  the  flight  path  may  be  written  In 
the  nomenclature  of  figure  1  as 


1  dV  -CnP^Ar  -CrnV^Agr 


(1) 


X  ^  +  (v2  .  i)cos  (P  »  +  (v2  -  Dcos  <p 


2W 


(2) 


The  symbols  are  defined  in  appendix  A. 

The  following  assumptions  are  made:  (l)  that  the  flight-path  angle 
relative  to  the  local  horizontal  is  always  small,  so  that  cos  cp  Is  ap¬ 
proximately  unity  and  sin  9  is  approximately  9,  and  (2)  that  sin  9 
Is  small  compared  with  the  nondlmenslonal  term  representing  the  aerody¬ 
namic  forces  in  equation  (l) . 


Equations  (l)  and  (2)  then  reduce  to 


1  ^ 
g  dt 


V  ^ 
g  dr 


1 


(3) 

(4a) 


where  G  is  defined  as 


CrpV^At  _  Opp^Agr 
2m  ^  2W 


(5) 
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G  represents  the  net  force  on  the  vehicle  in  terms  of  the  planet  g's, 
or  is  the  net  acceleration  vector  felt  by  the  pilot  in  planet  g's. 

For  future  reference,  equation  (4a),  with  the  aid  of  equation  (s) 
and  the  definition  of  (appendix  A),  may  also  be  written  as 


Also  from  figure  1,  the  component  of  lift  in  the  vertical  direction  is 
related  to  the  total  aerodynamic  lift  roll  angle  9  by 

=  ^L,a  9  (6) 

Vehicle  Aerodynamics 

The  vehicle  aerodynamics  (fig.  2)  are  calculated  using  the  equations 
and  assumptions  given  in  appendix  B.  The  vehicle  design  is  specified  by, 
or  characterized  by,  its  maximxan  lift-drag  ratio  In  general, 

for  a  given  lift  coefficient  these  curves  have  double  values  of  Cp,  Cp, 
and  a.  In  the  present  study  the  low  range  of  angle  of  attack  is  used. 

Flight  Paths 

Memy  possible  flight  paths  can  be  considered;  the  present  report 
considers  six  particular  flight  paths  and  a  generalized  flight  path.  The 
flight  paths  are  defined  as  follows: 

(1)  Constant  a:  0  =  0;  Cj^,  Cjj/Cp,  and  are  constant  (see 

fig.  2). 

(2)  Constant  h:  0  =  0;  p  is  constant. 

(3)  Constant  G:  0  =  0;  is  constant  (see  eq.  (s)). 

(4)  Constant  Re;  0  =  0;  pV  is  constant  (see  eq.  (17)). 

(5)  Constant  q:  0  =  0;  p  '  V  is  constant  (see  eq.  (lO)). 
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(6)  Modulated  roll;  6  is  varied  from  0°  to  90°  (or  180°)  in  a  man' 
ner  to  satisfy  equations  (4)  and  (6);  the  vehicle  angle  of  attack  is  de¬ 
fined  to  he  constant,  and  hence  Cp,  and  constant. 


Heat  Treuisfer 


T?he  heat  transfer  for  the  present  analysis  is  estimated  assuming  a 
laminar  boundary  layer.  References  8  and  9  present  a  comprehensive  re¬ 
view  of  heat-transfer  theories.  These  have  been  found  adequate  up  to 
Earth  circular  speeds  and  are  assumed  herein  to  be  applicable  also  at 
supercircular  speeds. 

Reference  10  gives  the  following  form  for  stagnation-point  heat 
transfer : 


15.5X10“%<2J/^pV^V^r  ^ 

4^2  ■  V 


(7a) 


where  J  =  0  for  the  two-dimensional  leading  edge  considered  in  the 
present  analysis  (j  =  1  for  a  three-dimensional  or  spherical  leading 
,  h^ 

edge) .  For  entry  from  high  velocities  «  1.0  and  with  equation 

“se 


(7a)  in  a  form  using  V, 


'is 


15 .5x10"3 


(7b) 


This  equation  may  be  modified  to  include  the  effect  of  sweep  and  nose 
radius  angle  by  the  relation  suggested  in  reference  8  (pp.  19  and  22): 


—  =  cos'^Aj,  cos  A  =  ^  -  sin^A  cos^c^'^^^  cos  A 


Combining  equations  (?)  eind  (8)  gives 

q,  =  15.5Xl0"^(gr)V2  P-  cos'^A 

Ry2 


cos  A 


(e) 


(9) 


The  theory  and  experiments  of  reference  11  indicate  that,  for  the  high 
velocities  of  the  present  study,  v  has  a  value  of  1 .5  for  sweep  angles 
from  0*^  to  60°  and  becomes  less  than  1.5  for  sweeps  greater  than  60°. 
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(The  theory  and  experiment  of  ref.  12  indicate  that  the  value  of  v  may 
be  as  low  as  unity.) 

The  symbols  defining  the  vehicle  geometry  and  various  flow  compo¬ 
nents  are  illustrated  in  figure  3.  In  calculating  the  heating  rates  and 
total  heat  inputs,  a  family  of  vehicles  was  considered  which  are  delta- 
wing  configurations  having  a  wing  sweepback  angle  A,  a  constant  thick¬ 
ness  d,  and  a  constant  nose  radius  Rjj.  The  configurations  have  no 

separate  fuselage,  and  the  volume  in  the  wing  is  considered  the  usable 
volume. 

The  following  analogously  simple  equation  for  the  heat  transfer  to 
the  windward  side  of  the  vehicle  is  developed  in  appendix  C: 


q„  =  15 .OXlO’^k^gr  Va  at  ^  =  10 . 2  kp  Va  aT  ( 10) 

c»  f  Ax  cos  6  °  f  Ax  cos  6  '  ' 


This  result  is  Independent  of  angle  of  attack  between  about  3°  and  30° 
for  a  constant  roll  angle  0. 


General  Form  of  Equations  for  Flight-Path  Characteristics 

In  addition  to  those  characteristics  involved  directly  with  the 
equations  of  motion  there  are  other  flight-path  characteristics  of  in¬ 
terest  and  these  are  given  in  general  form  here. 

Altitude .  -  With  the  assumption  of  an  exponential  atmosphere,  the 
flight-path  altitude  may  te  calculated  if  the  flight-path  density  is 
known; 


h  =  -  I  In  ^  (11) 

P  PSL 

For  Earth  the  values  =  2.35X10^  feet  and  Pgj^  =  0.0027  (slug)  (ft “2) 

are  used  as  suggested  in  reference  4.  Also  for  later  use  it  is  noted 
that 


dp  “  ■  3p 


(12) 


Radiation  equilibrium  temperature.  -  With  the  flight  path  and  ve- 
hicle  design  knovm,  equations  (9)  and  (lO)  give  the  heating  rate  per 
square  foot  q.  In  same  applications  the  radiation  equilibrium  tempera¬ 
ture,  determined  by  the  heat  rate,  is  of  interest;  it  is  given  by 


^w  “ 


1/4 


(13) 
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Free-stream  Reynolds  number.  -  With  the  flight -path  density  p 
luiovm  as  a  function  of  ?  and  the  vehicle  characteristic  length,  the 
free-stream  Reynolds  number  may  be  calculated  from 

Re  =  (14) 


By  use  of  equation  (5),  equation  (14)  may  alternatively  be  written  as 


Re 


(15) 


The  free-stream  viscosity  is  assumed  constant  and  Independent  of  altitude 
for  the  present  analysis,  p  =  0.37x10®  slugs  per  foot-second,  and  x  is 
taken  as  f  plus  the  root  chord  c  (see  fig.  3). 

Local  Reynolds  number.  -  Local  Reynolds  number  is  related  to  the 
free-stream  Reynolds  number  by  a  function  of  the  free-stream  Mach  number, 
the  wing-sweep  angle,  and  the  angle  of  attack.  The  method  of  calculation 
of  the  ratio  of  local  to  free-stream  Reynolds  number  Rej/Re  is  given 

in  appendix  D  and  assumes  T  =  1-2.  The  local  Reynolds  number  is  then 


The  value  of  Rej/Re  depends  on  the  value  selected  for  y,  and  more 

generally  on  the  real  gas  properties.  A  more  precise  determination  of 
the  local  Reynolds  number  than  that  of  appendix  D  is  of  little  value  at 
present,  because  the  treuisltion  Reynolds  number  is  not  accurately  known. 

Flight -path  angle.  -  As  mentioned  previously,  it  is  assumed  that 
the  flight-path  angles  are  small,  and  hence 

V 

sin  cp  ~  cp  -  (17) 


But 


V  dh  dh  dp  dV 
V  “  dt  ~  dp  dV  dT 


and  by  use  of  equations  (3),  (5),  and  (l2) 


9  = 


2P  j  dV 


(18) 


(19) 
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The  preceding  characteristics  are  point  functions  along  the  flight 
path.  The  succeeding  characteristics  result  from  an  integration  along 
the  flight  path. 


Time  of  flight.  -  The  time  of  flight  is  obtained  by  integrating  the 
time  increment  from  equation  (3)  as  follows; 


Range .  -  Using  equation  (3)  for  dr  and  approximating  cos  cp  by 
unity  give  the  range  traversed  during  the  deceleration: 


Total  heat  input.  -  One  of  the  important  criteria  of  merit  for  a 
heat -protection  system  is  its  weight  as  a  fraction  of  the  vehicle  gross 
weight.  If,  for  example,  the  heat-protection  system  is  of  the  ablation 
type  and  has  a  characteristic  heat-absorbing  capacity  C  in  Btu  per 
pound,  then  the  ablated  weight  may  be  expressed  as 


1  a 

w  ~  c  w 


(22) 


The  present  analysis  deals  with  the  parameter  Q/w  which  results  frcm 
a  double  Integration  (along  the  flight  path  and  over  the  vehicle) .  The 
Integration  along  the  flight  path  yields  the  total  heat  input  per  square 
foot,  which  is  given  by 


(23) 
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where  q  is  given  by  equations  (9)  and  (lO)  and  d-r  Is  again  given  by 
equation  (3) . 

The  integration  over  the  vehicle  has  the  general  form 


%  +  I 


q^Rjj  dA  dZ 


(24) 


This  integration  is  discussed  in  detail  in  appendix  E. 


Sample  Development  of  Flight-Path  Characteristics 
for  a  Constant  Reynolds  Number  Path 

The  case  of  "constant  free-stream  Reynolds  number"  is  considered 
here  in  detail.  The  velocity  ratio  V  is  taken  as  the  independent  vari¬ 
able.  The  free-stream  Reynolds  number  is  given  by  equation  (14)  or  (15),  * 
if  the  vehicle  is  specified  (i.e.,  f+c  and  w/A  are  constant  and  given). 
Then  for  Re  to  be  constant,  pV  must  be  consteuit.  Since  G  (eq.  (5)) 
and  all  the  terms  collected  on  the  left  side  of  the  equation 


are  defined  as  constants,  the  right  side  must  also  be  a  constant.  This 
is  the  relation  given  in  column  2  of  table  I.  The  first  column  identifies 
the  type  of  flight  path.  The  term  in  braces  is  a  constant  for  the  flight 
path  and  is  specified  if  the  path  Reynolds  number  is  specified. 


give  column  3  in  table  I: 


solved 

for 

1 

r 

2 

[oj 

^1 

A  gr 

V 

as  a  function  of  V  to 


(26) 


Because  {g/CrV}  is  a  consteuit  along  the  flight  path. 


(27) 
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which  is  the  result  given  in  column  4.  To  determine  the 
angle  cp,  the  derivative  dp/dV  is  required  in  equation 
be  obtained  from  equation  (26): 


flight -path 
(19) .  This  may 


(28) 


Combining  equations  (19)  and  (28)  leads  to  the  result  in  column  5: 


Cd 

[aJ 

1 

3r' 

V 

(29) 


With  cp  having  been  determined,  the  term 
equation  (29) 


—  ^  can  be  discussed.  From 
g  dT 


V  ^ 

g  dT  3r 


r  ^ 

G 


IPrV 


i.  1  ^ 

V  g  dT 


(30) 


where  the  term  - 

g  dT 

is  given  by  equation 
tion  (30)  becomes 


is  given  by  equation  (3),  and  G  in  that  equation 
(27)  (or  column  4).  With  these  substitutions  equa- 


(31) 


The  term  ^  ^  (column  6)  is  equal  to  ~%e  (column  19).  For  Cp/Cj^ 
and  Cp  constant,  an  assumption  required  later,  l^g  is  constant. 


Equation  (4a)  may  now  be  solved  for  the  lift  coefficient,  using  the 
result  of  equation  (3l)  and  equation  (27)  for  G,  to  arrive  at  the  re¬ 
sult  in  column  7 : 


(32) 


The  free-stream  Reynolds  number  (column  8)  is  given  by  equation  (14) 
and  is  a  constant  for  thl£  case.  More  generally  it  may  be  obtained  from 
the  relations  of  p  to  V  like  equation  (26)  or  column  3  of  table  I. 
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Obtaining  the  local  Reynolds  number  (column  9)  requires  knowledge 
of  the  wing  sweep  angle  and  the  angle  of  attack.  The  wing  sweep  is  part 
of  the  vehicle  design  specification.  The  angle  of  attack  is  determined 
by  figure  2  using  from  equation  (32)  and  the  vehicle  (lj/®)niax* 

The  heating  rate  to  the  leading  edge  q^  (eq.  (9))  now  depends  only 

on  the  parameter  p,  which  has  been  determined,  and  on  the  Independent 
variable  V.  This  result  is  shown  in  column  10.  The  heating  rate  to 
the  windward  side  q^  (eq.  (lO))  depends  on  L/A,  which  may  be  determined 

from  equation  (4b)  and  the  result  of  equation  (31)  or  column  6  as  fol¬ 
lows: 

I  =  I  (1  -  V2  -  kRe)  (33) 

Making  this  substitution  in  equation  (lO)  gives  the  result  listed  in 
column  11  of  table  I. 

The  point  functions  have  now  all  been  determined,  and  it  remains 
to  determine  the  Integrated  functions.  The  time  of  flight  may  be  ob¬ 
tained  from  equation  (20)  and  the  equation  for  0  from  column  4  (or 
eq.  (27)),  where  the  term  in  braces  is  a  constant,  so  that 


The  variation  of  ^  Just  Cp)  as  a  function  of  V  re¬ 

quired  in  equation  (34)  can  be  obtained  from  the  curves  of  figure  2  and 
equation  (32)  for  Cj^.  Hence,  equation  (34)  can  be  integrated  graphi¬ 
cally  or  by  numerical  methods  for  the  most  general  cases.  The  present 
report,  however,  deals  only  with  the  more  simple  cases  where  Cp/Cj^  and 

Cp  can  be  approximated  by  constant  values  and  closed-form  solutions  can 

be  thus  obtained.  This  assumption  is  different  from  the  assumptions  for 
the  case  of  constant  a  for  which  C^/Oj^  in  addition  to  Cp/Cj^  and 

Cp  is  constant.  With  the  assumptions  of  Cp/Cp  and  constant. 
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which  are  noted  in  column  12,  equation  (34)  becomes 


which  integrates  to  the  result  in  column  13. 

By  similar  substitutions  and  assumptions,  equation  (2l)  for  this 
range  becomes 


which  integrates  to  the  result  in  column  14. 


Again  with  the  similar  substitutions  and  approximations,  and  by 
use  of  the  result  in  column  10  or  equation  (9)  for  q^,  equation  (23) 

for  q2  becomes 


which  integrates  to  the  result  in  column  15. 

When  column  11  or  equation  (lO)  with  0=0  is  used  for  q^  in 
equation  (23),  the  integral  equation  becomes 


q 


w 
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which  integrates  to  the  result  in  column  16. 

The  integration  over  the  vehicle  to  obtain  q/W  indicated  by  equa¬ 
tion  (24)  is  carried  out  in  detail  in  appendix  E  and  gives  the  results 
shown  in  columns  17  and  18.  Alternative  forms  for  the  denominator  of 
columns  17  to  18  are  shown  in  the  bottom  row. 


E-lOOl 


The  results  for  all  the  other  flight  paths  except  constant  a  can 
be  derived  in  a  similar  manner-  The  case  of  constant  angle  of  attack  is 
discussed  in  appendix  F. 


RESULTS  AND  DISCUSSION 

The  results  of  the  analysis  for  all  the  flight  paths  studied  are 
summarized  in  tables  I  and  II. 


Tables  of  Equations 

The  first  column  of  table  I  gives  the  name  of  the  flight  path.  The 
term  in  braces  in  the  second  and  subsequent  columns  is  constant  over  the 
flight  path.  The  equations  with  the  exception  of  those  in  columns  10, 

11,  and  15  to  18  apply  to  any  planet.  The  results  in  the  aforementioned 
columns  are  given  in  table  II  in  a  form  which  is  applicable  to  any  planet 
that  has  an  atmosphere  with  heat-transfer  characteristics  similar  to 
Earth.  The  equations  for  G  are  in  terms  of  planet  g's,  eind  V  is  the 
flight  velocity  divided  by  the  planet  circular  velocity. 

Three  heating  parameters  are  presented: 

(1)  Heating  rate  q  (in  Btu/(sq  ft)(sec)).  This  is  the  parameter 
of  interest,  for  example,  for  radiation  cooling. 

(2)  Heat  input  per  square  foot  of  vehicle  q,  which  depends  on  the 
location  with  respect  to  the  stagnation  point  of  the  area  being  con¬ 
sidered.  (The  stagnation  point  is  not  a  fixed  geometric  point  except 
for  the  cases  of  constant  a  and  modulated  roll.)  The  peurameter  q  re¬ 
sults  from  the  integration  over  the  flight  path  of  the  heat  input  to  a 
given  location. 

(3)  Heat  input  Q  =  (Qj  +  f  which  is  related  to  the  ablated  weight 

and  involves  an  integration  over  the  vehicle  surface  as  well  as  over  the 
flight  path. 

The  heating  parameters  are  presented  for  two  regions  on  the  airplane : 
(1)  the  region  in  the  vicinity  of  the  leading  edge  or  stagnation  point, 
where  the  heating  is  a  function  of  nose  radius  Rj^,  nose  radius  coordi¬ 
nate  angle  A,  and  sweep  angle  A,  and  (2)  the  windward  side  of  the  ve¬ 
hicle,  where  the  heat  input  is  a  function  of  the  distance  x  from  the 
stagnation  point.  The  heat  input  to  base  and  lee  regions  is  not  con¬ 
sidered. 
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The  term  k2  (column  17)  accounts  for  the  reduction  with  angle  of 

attack  of  both  the  effective  sweepback  angle  and  the  area  of  the  leading 
edge  subject  to  heating.  The  value  of  k2  is  approximately  unity.  The 

term  k^,  (columns  11,  16,  and  18)  accounts  for  the  relative  lift  generated 
by  the  windward  side  and  the  leading  edge.  The  term  k,^  accounts  for 
the  length  of  boundary  layer  ahead  of  the  windward  side.  The  term  k 
•with  a  subscript  to  designate  the  flight  path  (e.g.,  kj^g)  accounts  for 

the  amount  that  the  flight-path  curvature  differs  from  the  planet  surface 
curvature.  Both  k^  and  k^  approach  unity  for  slender  (or  high 

(l/d)^^^)  configurations  and  have  a  value  of  about  0.5  for  the 
(L/D)jjjg^  =1.0  configuration  of  this  report. 

In  table  I  the  heading  "Modulated  angle  of  attack"  indicates  the 
flight  paths  for  which  the  lift  vector  is  constrained  to  the  vertical 
direction.  For  all  but  the  constant  a  case,  the  assumption  of  con¬ 
stant  Cd/C|2  restricts  the  angle  of  attack  to  the  vicinity  of  0°  or  90°, 

as  may  be  seen  from  figure  2.  The  present  results  are  for  a  near  0°. 

"Modulated  roll"  in  table  I  designates  a  specific  case  where  angle 
of  attack  with  respect  to  the  free-stream  direction  is  constant  (hence, 

Cp,  Cj)/Cp,  and  constant)  and  the  variation  of  the  roll  ang].e 

6  with  Jf  is  given  by  the  equation  in  column  1^  This  relation  between 

d  and  V  gives  a  constant  G  and  0  =  0  at 


Comparison  with  Other  Analyses 

To  obtain  closed-form  solutions,  a  number  of  simplifying  assumptions 
have  been  made.  As  an  indication  of  the  accuracy  of  the  present  analysis 
in  the  stagnation  region,  a  comparison  is  made  for  the  constant-angle- 
of-attack  path  with  the  results  from  machine  integrations  of  the  equa¬ 
tions  of  motion  and  heating  presented  in  reference  4.  The  procedure  for 
reducing  the  results  of  the  present  report  to  the  form  of  reference  4  is 
given  in  appendix  G.  Figure  4  presents  the  horizontal  component  of  the 
G  load  Oh  and  the  heating-rate  parameter  q  as  a  function  of  velocity 

ratio  V  for  an  L/D  of  0.1  and  1.0  and  an  initial  velocity  ratio  V 
of  0_^99.  The  heat-input  parameter  q  for  decelerating  from  V  =  0j_99 
to  V  =  0  is  given  as  a  function  of  L/D.  In  the  calculation  of  £, 
the  function  F  of  table  I  was  evaluated  at  the  maximum  value  of  V. 

The  present  analysis  generally  gives  the  greater  heat  input.  In  all 
cases  the  agreement  between  the  two  analyses  is  considered  satisfactory. 
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To  indicate  the  accuracy  of  the  present  analysis  in  estimating  the 
heat  input  to  the  windward  side  of  the  vehicle,  a  comparison  is  made  with 
the  results  of  reference  13  for  the  constant-angle-of-attack  path.  The 
results  of  reference  13  have  been  extrapolated  to  the  condition  of  zero 
wall  temperature  and  apply  for  the  case  of  P  =  1  in  table  I .  Figure  5 
compares  the  heating-rate  parameter  4(x/(w/A) and  the  heat  input 

•c' 

q  dx 

'0 

as  a  function  of  velocity  ratio  V 


per  unit  span  parameter 


/ 

L/c'wy 

dU  / 


•wW^ 


for  an  initial  velocity  V  of  1.0.  Also  shown  is  the  variation  of  heat¬ 
ing  rate  with  angle  of  attack  for  constant  lift.  In  general,  the  agree¬ 
ment  between  the  two  aneJ-yses  is  satisfactory.  Compared  with  reference 
13,  the  present  analysis  indicates  Increasingly  higher  heat  inyuts  with 
increasing  angle  of  attack.  The  maximum  angle  of  attack  to  which  the 
present  analysis  is  applicable  depends  on  the  degree  of  approximation 
considered  acceptable.  The  present  analysis  is  about  15  percent  higher 
at  a  =  30°. 


In  general,  the  accxiracy  of  the  present  analysis  is  adequate  for 
preliminary  studies. 


Seme  General  Considerations 

For  atmospheric  bredting  it  is  desired  to  determine  a  combination 
of  vehicle  geometry  and  flight  path  which  will  result  in  a  low  heat  in¬ 
put  to  the  vehicle.  Generally  there  is  a  substantial  increase  in  con¬ 
vective  heat  input  to  the  vehicle  if  the  boundary  layer  on  the  vehicle 
is  turbulent  rather  than  laminar.  Thus  it  is  desirable  to  restrict  the 
flight  path  to  conditions  for  which  laminar  boundary  layer  can  be  main¬ 
tained.  An  indication  of  the  possible  nature  of  the  boundary  layer  is 
given  by  the  magnitude  of  the  Reynolds  number.  Two  Reynolds  nianbers  ceui 
be  considered,  the  free-stream  value,  and  a  local  value  at  the  edge  of 
the  boundary  layer.  The  free-stream  Reynolds  nimiber  depends  on  the 
vehicle  design  parameters  W/A  and  c  and  on  the  flight-path  parameter 
G  (see  eq.  (17)  or  (18)).  The  local  Reynolds  number  is  further  influ¬ 
enced  by  the  vehicle  design  factors  of  leading-edge  bluntness  and  sweep 
and  the  operating  parameter  of  angle  of  attack.  Both  local  and  free- 
stream  Reynolds  number  are  presented  in  the  section  Point  Functions. 

The  local  Reynolds  number  is  interrelated  with  the  problem  of  hot- 
gas  radiation.  At  entry  velocity  ratios  V  greater  than  1.0,  the  heat 
input  from  hot-gas  radiation  can  equal  or  become  larger  than  the  con¬ 
vective  heat  input  to  an  unswept  stagnation  region.  The  lowest  heat  in¬ 
put  to  the  vehicle  will  usually  exist  if  hot-gas  radiation  is  avoided 
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where  possible .  This  requires  that  the  vehicle  surface  angles  be  small 
with  respect  to  the  free-stream  direction  to  avoid  high  air  static  tem¬ 
peratures.  In  general,  unsweeping  the  wing  leading  edge  decreases  the 
local  Reynolds  number  but  tends  to  aggravate  the  problem  of  hot-gas  radia¬ 
tion.  The  choice  of  wing  sweep  angle  is,  therefore,  a  compranise.  The 
calculations  for  the  present  study  were  made  for  a  sweep  angle  of  60°, 

Also  to  avoid  hot -gas  radiation  the  low-angle -of -attack  rather  than  high- 
angle-of -attack  region  of  aerodynamic  operation  was  chosen.  The  choice 
of  the  low-angle-of-attack  range  and  the  60°  sweep  is  assumed  to  reduce 
the  heat  input  from  hot-gas  radiation  to  a  low  value  compared  with  the 
convective  heat  input.  Only  the  convective  heat  input  is  considered 
subsequently  in  this  report. 


Point  Functions 

The  evaluation  of  the  various  flight  paths  for  decelerating  in  the 
Earth  atmosphere  is  discussed  with  the  aid  of  figures  6  to  11  in  which 
is  plotted  the  variation  of  the  point  functions  of  table  I  with  V.  A 
^econd  abscissa  shows  elapsed  time  in  minutes  with  time  equed.  to  zero  at 
V  =  2.0.  To  arrive  at  numerical  results,  a  vehicle  design  defined  by 
Pj  =  20  pounds  per  cubic  foot,  W  =  10,000  pounds,  A  =  60°,  and 

(L/D)jugjj  =1.0  was  assumed.  (Higher  values  of  (L/d)jj,^  result  in 

larger  heat  inputs.)  For  a  vehicle  like  that  of  figure  3  and  the  above 
numerical  values,  the  following  vehicle  characteristics  may  be  determined 
by  the  equations  given  in  appendix  E;  W/a  =  83  pounds  per  square  foot, 

Rjj  =  2.06  feet,  and  c  =  14..‘i  feet.  A  length  of  18.8  feet  was  used  in 

the  calculation  of  the  Reynolds  numbers.  Vehicle  operation  is  in  the 
low-angle-of-attack  reuige,  and  the  assumptions  are  made  that  Cjj 

and  Cp/Cp  =  1.0,  although  for  some  cases  more  precise  approaches  to  the 

calculation  of  the  point  functions  are  possible.  All  the  flight  paths 
have  a  common  point,  at  V  =  2.0,  G  =  10  g's,  and  9=0.  This  results 
in  approximately  the  same  values  of  a  and  h,  about  5.3°  and  162,000 
feet,  respectively.  When  appropriate,  the  following  physical  limits  are 
noted  on  the  curves:  |CL^niaxl  “  ®max  ~  S's,  hjjin  = 

4  =  190  Btu/(ft"^) (sec) ,  or  t^  =  4460°  R  for  an  emissivity  of  unity. 

Flight -path  angle.  -  The  fact  that  for  all  cases  the  flight-path 
angle  9  is  small  except  in  the  vicinity  of  V  =  0  justifies  the  assump¬ 
tion  of  small  angles.  For  all  cases  except  constant  a  (fig.  6),  9  is 
also  continuous.  For  the_constant  a  case,  9  changes  discontinuously 
from  plus  to  minus  near  V  =  1.0  and  the  altitude  becomes  indefinite. 

This  situation  j^s  physically  approached  by  a  shallow  skip  out  of  the 
atmosphere  at  V  =  1.0.  In  spite  of  this  indefiniteness  in  the  altitude, 
the  heat  input  to  the  vehicle  remains  finite  as  is  shown  later  (see  also 
columns  17  and  18  of  table  l) . 
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Altitude  and  lift  coefficients.  -  For  all  cases  except  constant  h 
(fig.  7)  the  altitude  becomes  less  than  zero  for  V  approaching  zero. 

On  the  constant  h  path  excessive  lift-coefficient  requirements  are  en¬ 
countered  at  V‘s  approaching  zero.  The  constant_  Re  path  (fig.  8) 
also  requires  excessive  lift  coefficients  at  low  Vs.  In  a  practical 
flight  these  and  other  low-speed  anctnalies  are  usually  precluded  by 
specifying  some  terminal  condition  of  velocity,  lift  coefficient,  or 
altitude .  The  actual  flight  path  may  deviate  from  the  specified  path  in 
the  low-speed  range  with  little  effect  on  the  total  heat  input  calculated 
from  an  initially  supercircular  velocity,  as  is  shown  later. 

The  lift  coefficient  is  negative;  that  is,  lift  force  is  directed 
toward  the  Earth  for  flight  velocities  greater  than  circular.  This 
occurs  fundamentally  because  the  centrifugal  force  tending  to  "throw" 
the  vehicle  out  of  the  atmosphere  is  greater  than  the  gravity  force  and 
the  difference  in  these  forces  must  be  counteracted  by  the  lift.  For 
all  cases  except  constant  G  (fig.  9),  constant  q  (fig.  10),  and  modu¬ 
lated  roll  (fig.  ll),  the  lift  coefficient  is  positive  (directed  away 
from  the  Eeurth)  for  subcircular  velocities,  because  at  subcircular  speed 
the  centrifugal  force  is  less  them  the  gravity  force.  However,  to  satisfy 
the  assigned  veilues  of  q  and  Q  for  the  cases  of  constant  q,  constemt 
G,  and  modulated  roll,  negative  lift  coefficients  are  required  at  sub- 
circular  speeds  as  well  as  at  supercircular  speeds  (because  of  the  term 
(V/g)(d<p/dt)  in  eq.  (4)). 

G  load.  -  For  the  constant  G  path  (G  =  10  g's),  the  abscissa  of 
figure  9  shows  that  the  time  of  deceleration  in  the  region  of  strong 
heating,  that  is,  from  V  =  2.0  to  V~0.5,  is  about  2  minutes.  This 
combination  of  g's  and  time  is  about  the  presently  accepted  limit  of 
human  tolerance  if  useful  functions  are  to  be  performed  by  the  mem.  For 
this  reason  G^iax  “  used  in  the  present  examples. 

For  all  cases  except  constant  4  (fig*  IG)  the  G  load  is  at  or 
below  the  initial  G  value  at  V  =_2.0.  For  the  case  of  constant  q 
the  G  increases  with  decreasing  V  to  maintain  q  constant.  At  a 
V  higher  than  2  the  Cj^  becomes  excessive.  These  same  characteristics 

exist  if  other  values  of  constant  q  are  specified.  Thus  it  is  con¬ 
cluded  that  the  constant  q  path  does  not  have  application  over  a  wide 
range  of  V . 

Angle  of  attack.  -  The  angle  of  attack  is  shown  as  positive.  For 
the  flight  paths  where  the  lift  coefficient  varies  from  minus  to  plus 
along  the  flight  path,  one  has  the  choice  of  permitting  the  sign  of  the 
angle  of  attack  to  change  or  of  rolling  the  vehicle  180°.  In  the  latter 
instance  the  same  side  of  the  vehicle  remains  the  windward  side.  The 
vehicle  angle  of  attack  is  for  all  cases  equal  to  or  less  than  the  angle 
of  attack  at  V  =  2.0,  where  a  is  about  5.3°  (except  for  the  very  low 
Vs).  From  figure  2  it  can  be  seen  that  Cj^,  Cp/Cj^,  and  vary  less 
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than  ±5  percent  from  a  mean  value  over  the  angle  of  attack  reuige  0°  to 
5.3°.  More  generally,  for  =  10  g's  and  (0^/0^)^  =  (l  -  V^)/G^ 

^eq .  (4a)  with  ^  =  o)  >  Cjj,  for  example,  varies  less  than  ±7  percent 
from  a  mean  value  up  to  =  2.25  for  (l/d)j^qjj  =  1  and  2,  the  values 
of  (l'/^)niax  examined.  An  increase  in  Gj_  will  increase  Vj^,  for  which 

the  variation  of  the  coefficients  is  small.  Hence,  for  many  problems  of 
interest  it  is  a  good  approximation  to  use  constant  values  for  Cj), 

Cjj/Cj^,  and  as  is  required  to  arrive  at  some  of  the  integrated  func¬ 

tions. 

IVee-stream  Reynolds  nxmiber.  -  The  free-stream  Reynolds  number  Re 
is  slLso  shown  in  figures  6  to  11.  The  case  of  constant  h  shows  a 
linearly  decreasing  Reynolds  number  with  decreasing  V.  The  case  of  con¬ 
stant  G  for  G  *  10  g's  shows  increasing  Reynolds  number  with  decreasing 

V  and  values  much  larger  than  for  the  other  cases  in  the  lower  V  range . 
The  various  flight  paths  Illustrated  thus  offer  a  choice  of  increasing, 
constant,  or  decreasing  Reynolds  number  with  decreasing  V. 

Local  Reynolds  number.  -  The  transition  from  laminar  to  turbulent 
flow  in  the  boundary  layer  is  more  significantly  related  to  the  loced 
Reynolds  number  than  to  the  free-stream  Reynolds  number.  The  variation 
of  the  local  Reynolds  number  Rej  with  V  is  shown  for  the  condition 

that  the  wing  leading  edge  is  blunt  and  swept  60°.  All  the  flight  paths 
except  modulated  roll  (fig.  11)  show  a  dip  in  Re^  in  the  vicinity  of 

V  =  1  followed  by  a  rise  in  Re^  for  V  <1.0.  For  the  flight  paths  of 

constant  G  and  constant  Re,  this  dip  is  associated  with  a  lower  angle 
of  attack  in  the  vicinity  of  V  =  1.0.  For  the  constant  a  case  it  is 
associated  with  a  low  air  density  (see  fig.  6).  In  the  vicinity  of 

V  =  1.0,  the  local  Reynolds  number  for  both  the  constant  a  (a  ~  5.3°) 
and  constant  G  (G  =  10)  cases  is  svibstemtially  lower  them  for  the 
modulated  roll  case  (G  =  10  g's  and  a  =  5.3°).  This  illustrates  e.  dis¬ 
advantage  of  using  modulated  roll.  The  high  local  Reynolds  number  for 
modulated  roll  is  associated  with  the  £ombination  of  the  high  air  den¬ 
sity  and  the  5.3°  angle  of  attack  at  V  =  1.0.  The  effect  of  angle  of 
atteu:k  on  the  ratio  Re^/Re  is  discussed  in  appendix  D.  In  general, 

for  the  flight  paths  considered,  it  appears  less  difficult  to  achieve 
low  local  Reynolds  number  at  supercircular  than  at  subcircular  speeds. 

The  constant  a,  constant  h,  constant  Re,  and  constant  G  flight 
paths  show  rather  small  variation  in  Rej  with  V  and  little  variation 

among  the  several  flight  paths  over  the  speed  range  from  V  =  3.0  to 
0.8.  In  this  speed  range  the  general  level  of  the  local  Reynolds  number, 
0.7X10^,  suggests  an  extensive  run  of  laminar  boundary  layer  on  the  ve¬ 
hicle  is  possible. 
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Heating  rate.  -  The  values  of  the  heating  rate  q  shewn  are  for  a 
vehicle  density  Pf  =  20  (lb)(ft“3)  and  apply  to  the  swept  leading  edge 

(A  =  60°)  (note  that,  in  column  10  of  table  I,  W/RjjA  =  2pf ) .  For  eLLl 

cases  except  constant  4»  values  of  4  Increase  with  increasing  V . 

The  temperature  limits  of  presently  known  materials  (e.g.,  ceramics)  are 
about  4460°  R.  By  equation  (16)  this  corresponds  to 

q  =  190  (Btu) (ft“2) (sec"^) .  For  this  value  the  best  case,  constant  a, 
shows  that  entry  velocity  ratios  V  up  to  1.3  may  be  possible  with 
radiation  cooling,  for  the  present  vehicle  and  flight  path.  (The  flight 
path  and  vehicle  design  shown  do  not  yield  the  lowest  possible  radiation 
y  temperatures.)  For  entiy  velocity  ratios  greater  than  approximately  1.4 

radiation  cooling  is  probably  not  possible,  and  hence  there  is  interest 
in  the  total  heat  input  to  the  vehicle. 

Maximum  entry  velocity.  -  Fran  equation  4(a)  for  =  10  g's  and 
the  maximum  values  of  shown  in  figure  2,  it  can  be  determined 

that  entry  velocities  of  2.8  times  circular  are  aerodynamically  possible 
with  a  vehicle  with  an  (L/D)ijjax  Still  higher  entry  velocities 

are  possible  with  a  vehicle  with  a  higher  (L/D)niax' 

Summary  of  point  functions.  -  The  more  important  conclusions  from 
the  discussion  of  the  point  functions  may  be  summarized  as  follows: 

(1)  All  the  flight  paths  considered  except  constant  4  or  ap- 

*  proximate  physically  possible  flight  paths. 

(2)  For  =  10  g’s  and  (l/d)^^^^  =  1.0  or  2.0  the  assumption  of 
Cp,  Cp,  and  Cjj/Cp  constant  is  a  valid  approximation  up  to  an  entry 
velocity  ratio  V  of  about  2.25. 

(3)  With  the  vehicle  designs  and  flight  paths  assumed,  radiation 
cooling  appears  Impossible  for  entries  at  supercircular  speeds  corre¬ 
sponding  to  V  >  1.4.  Hence,  the  total  heat  input  into  the  vehicle  is 
important . 

(4)  The  local  Reynolds  number  for  many  of  the  flight  paths  of  in¬ 
terest  is  about  0.7x10®.  This  value  suggests  the  possibility  of  a 
laminar  boundary  layer. 


Toted  Heat  Input 

For  the  assumptions  of  lamineu:  boundary  layer  and  hot-gas  radiation 
negligible  because  of  low  vehicle  surface  angles  and  with  heat  input  to 
base  and  lee  regions  neglected,  the  total  heat  input  to  the  vehicle  is 
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given  ty  the  sum  of  columns  17  and  10  of  table  I  or  II.  The  total  heat 
input  depends  on  many  factors,  and  these  may  be  categorized  as  follows: 

(1)  Factors  that  depend  on  the  flight  plan:_  terms  in  braces  or  in 
column  2  that  define  the  variation  of  G  with  V  and  the  terms  in 
brackets  in  columns  15  and  16  that  depend  on  the  initial  and  final  veloc¬ 
ity  of  the  deceleration  flight  path. 

(2)  Factors  that  depend  on  the  vehicle  design:  vehicle  density, 
gross  weight,  W;  wing  loading,  W/A;  and  wing  sweep  angle,  A. 

(3)  Factors  that  depend  on  both  vehicle  design  and  flight  plan:  Cjj 
and  L/D  or  Cp/Cp  and  Cj^.  These  depend  on  design  through  the  con¬ 
figuration  maximum  lift-drag  ratio  and  on  the  flight  plan  through  the 
operating  lift-drag  ratio  or  angle  of  attack. 

The  heat  input  along  families  of  flight  plans  related  in  two  ways 
to  the  curves  of  figures  6  to  11  will  be  considered.  For  the  first  family 
the  flight  plans  start  at  various  points  along  the  curves  of  figures  6  to 
11  and  the  flight  path  then  follows  these  curves .  The  second  family  is 
defined  by  the  condition  that,  at  the  initial  velocity  the  initial 

G  load  Gj^  is  10  g's. 

Figure  12  considers  the  first  family  of  flight  paths,  which  have  in 
common  a  net  acceleration  of  10  g's  at  V  =  2.0,  and  presents  the  heat 
to  the  windward  side  of  the  vehicle  ^  and  the  total  heat  input  ^ 

(the  sum  of  the  heat  inputs  to  the  windward  side  and  to  the  leading-edge 
region)  as  a  fiinction  of  Though  not  shown  by  the  curves,  from  the 

equation  of  tables  I  and  II,  the  heat  input  Q  equals  zero  for  =  0. 

Each  point  in  figure  12  can  be  interpreted  as  giving  the  heat  input 
for  a  complete  "deceleration  flight  path"  having  a  terminal  velocity  of 
zero  and  the  initial  velocity  and  initial  G  given  by  figures  6  to  11. 

The  history  of  the  flight  is  also  given  by  figures  6  to  11.  For  the 
supercircular  speed  range  the  flight  paths  in  order  of  increasing  total 
heat  input  for  a  given  Vj^  rank  as  follows:  constant  G,  modulated 
roll,  constant  Re,  constant  h,  and  constant  a.  At  =  2,  the  con¬ 
stant  a  case  has  about  twice  the  heat  input  of  the  constant  G  case. 

The  curves  of  figure  12  can  also  b£  interpreted  to  show  the  distri¬ 
bution  of  heat  input  as  a  function  of  V.  For  example,  for  a  flight 
starting  at  any  initial  velocity,  the  contribution  to  the  heat  input  in 
decelerating  from  V  =  0.5  to  0  is  given  by  the  Q  at  V  =  0.5.  With 
this  interpretation  for  all  the  flight  plans  the  heat  input  for  decelerat¬ 
ing  from  V  =  0.5  to  0  is  generally  an  order  of  magnitude  less  than  the 
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total  heat  input  for  decelerating  from  any  super circular  speed.  A  precise 
determination  of  the  total  heat  input  at  the  low  speeds  is  not  essentied 
to  an  accurate  estimate  of  the  total  heat  input  in  decelerating  from 
supercircular  speeds.  This  implies  that  neglecting  the  low-speed  anom- 
edies  previously  discussed  in  regard  to  the  point  functions  has  little 
effect  on  the  accuracy  of  calculating  the  total  heat  input  for  the  de¬ 
celeration  from  supercircular  speeds,  and  therefore  these  anomalies  have 
been  neglected  in  the  present  calculations.  For  most  flight  paths  the 
total  heat  input  to  the  windward  side  of  the  vehicle  is  about  15  percent 
of  the  total  heat  input,  for  a  vehicle  with  an  of  1.0. 

The  second  family  of  deceleration  flight  paths  considered  is  shown 
in  figure  13.  Each  point  in  the  figure  represents  a  complete  decelera¬ 
tion  flight  path  having  the  given  by  the  abscissa  and  •=  10  g's. 

The  values  for  =  2.0  for  all  the  flight  plans  and  for  the  complete 

curve  for  constant  G  are  the  same  for  figures  12  and  13.  In  the  super¬ 
circular  speed  range  the  fl'ght  paths  generally  rank  in  the  same  order 
as  for  the  preceding  family  of  curves.  The  equations  for  the  constant  a 
path  are  of  questj^onable  accuracy  for  the  very  low  L/D's  (l/D  <  0.03) 
ttat  restjit  for  Vj^  close  to  unity,  and  hence  results  are  not  shown  at 

V  less  than  1.2. 

The  relation  of  the  ablated  weight  to  the  total  heat  input  is  given 
by  equation  (26).  An  ablation  system  based  on  the  extrapolation  of  the 
data  of  references  14  and  15  to  the  conditions  at  supercircular  speeds 
and  a  system  utilizing  the  boiling  of  lithium  (ref.  16)  both  indicate 
the  possibility  of  an  effective  heat-absorbing  capacity  of  about  10,000 
Btu  per  pound  of  ablated  weight.  Using  this  value  for  0  in  equation 
(26)  for  =  2.0  and  G^^  =  10  g's  gives  values  of  the  ablated  weight 

as  a  fraction  of  the  vehicle  gross  weight  of  0.032  and  0.085  for  the 
flight  paths  of  constant  G  and  constant  a,  respectively.  Halving  the 
value  of  C  would  double  the  ablated  weight.  From  ;the  consideration  of 
the  ablated  weight  it  is  concluded  that  entries  at  V  =  2.0  and  greater 
are  feasible . 

The  distribution  of  heat  between  the  windward  side  and  the  leading 
edge  of  the  vehj  cle  depends  on  the  G  load,  as  illustrated  in  figure  14 
for  a  family  of  constant  a  flight  paths  with  Vj^  =  2.0  and  for  a  ve¬ 
hicle  with  an  (L/D)jjjgjj  of  1.0.  As  the  initial  G  load  decreases,  the 

vehicle  angle  of  attack  and  initial  flight  altitude  must  increase  to 
satisfy  equations  (4a)  and  (5).  Because  of  the  assumption  of  constant 
and  Cp/Cp  the  other  flight  plans  cannot  be  considered  over  this 

wide  a  range  of  G  or  a.  The  curves  show  that  the  heat  input  to  both 
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the  leading-edge  region  and  the  windward  side  of  the  vehicle  Increases 
with  decreasing  G  load.  Reducing  the  initial  Gj^  from  10  to  4.25  g's 

about  doubles  the  total  heat  input.  The  veilue  Gj^  =  4.25  g's  is  the 

minimum  at  which  the  constant  a  path  can  be  initiated.  Also  the  pro¬ 
portion  of  the  heat  input  to  the  windward  side  of  the  vehicle  increases 
fran  about  19  percent  at  Gj^  =  10  g's  to  47  percent  at  Gj^  =  4.25  g's. 

The  increased  angle  of  attack  accounts  for  the  increased  proportion  of 
heat  to  the  windward  side  of  the  vehicle. 

While  the  heat  input  for  the  constant  a,  G^  =  4.25  g  path  is  2.2 
times  that  for  the  constant  a,  Gj^  =  10  g  path,  it  is  4.9  times  larger 

than  that  for  the  constant  G  =  10  g  path.  This  indicates  the  importance 
of  both  the  G  level  euid  the  type  of  deceleration  flight  path  to  the 
vehicle  heat  input . 

Figure  14  also  shows  the  Reynolds  numbers  at  =  2.0.  Decreasing 

the  G  load  decreases  the  free-stream  Reynolds  number  almost  directly, 
as  would  be  expected  from  equation  (18)  (for  constant  Or).  Increasing 

angle  of  attack,  however,  increases  the  ratio  of  local  to  free-stream 
Reynolds  number  (see  appendix  D),  so  that  decreasing  Gj^  has  essentially 

no  effect  on  the  local  Reynolds  number. 


Control  of  Landing  Location 

One  of  the  practical  requirements  that  may  be  placed  on  an  entry 
vehicle  is  that  it  have  the  capability  of  landing  at  a  predesignated 
location.  The  problem  of  achieving  a  low  total  heat  input  on  a  flight 
path  that  will  also  yield  control  of  the  landing  location  can  be  studied 
to  a  first-order  approximation  with  the  present  solutions. 

To  land  at  a  predesignated  location  requires  control  of  motion  both 
down  and  cross  range  with  respect  to  the  Earth's  surface.  By  column  14 
of  table  I  the  down  range  for  the  constant  a  case  is  indefinitely  large. 
For  this  case  the  analysis  implies  that  at  V  =  1.0  the  vehicle  glides 
up  out  of  the  atmosphere  into  an  orbit  of  indefinite  duration.  (The  de¬ 
tails  of  this  maneuver  are  not  defined  by  the  equation.  If  the  vehicle 
exits  from  the  atmosphere  at  a  small  angle  as  indicated  by  fig.  6,  a 
small  rocket  may  be  required  to  "circularize"  the  orbit,  and  then  a  second 
small  rocket  may  be  required  to  initiate  reentry.)  Cross  range  movement 
may  be  obtained  by  waiting  for  the  Earth  to  rotate  through  the  plane  of 
the  orbit  If,  for  example,  the  vehicle  had  glided  into  a  polar  circular 
orbit,  all  of  Earth's  surface  passes  through  the  plane  of  the  vehicle  in 
1/2  revolution  of  Earth  or  in  12  hours.  This  time  corresponds  to  about 
8  revolutions  of  the  vehicle  about  Earth.  A  down-range  capability  of  8 


TAHT- 


E-lOOl 


23 


revolutions  thus  yields  a  control  of  cross  range.  With  the  eterodynamic 
cross-range  capability  such  that  half  the  distance  between  the  traces  on 
the  Earth's  surface  of  two  successive  orbits  can  be  traversed  during  the 
deceleration  from  the  circular  orbit,  the  vehicle  can  potentially  land 
at  any  predesignated  location. 

The  present  discussion  deals  with  the  down  range  and  total  heat  in¬ 
put  for  some  possible  flight  plans  that_include  control  of  landing  loca¬ 
tion  by  control  of  the  flight  through  V  =  1.0.  The  flight  plans  con- 
spidered  are  illustrated  in  figure  15  in  terms  of  the  G  history  against 
V .  The  flight  paths  are  composites  of  the  constant  G  and  constant  a 
flight  paths,  and  all  begin  at  V  =  2.0  with  an  initieLi  deceleration  of 
G  =  10  g's.  Figure  15(a)  shows  the  constant  a  case  modified  by  a  con¬ 
stant  G  segment  through  V  =  1.0.  The  magnitude  of  the  deceleration 
during  this  constant  G  segment  is  designated  because  it  is  the 

minimum  deceleration  that  occurs  during  the  deceleration.  Then  Gmin 
can  be  considered  an  independent  variable.  For  Gj^j^  =  0  the  general 
path  reduces  to  the  constant  a  case  and  for  Gmin  =  10  g's  to  the  con¬ 
stant  G  case . 

The  variation  of  Q  and  range  with  is  presented  in  figure 

16.  The  minimum  Q  and  range  occur  for  Gmin  “  values  are 

390  Btu  per  pound  of  gross  weight  and  0.032  of  the  Eeurth  circumference, 
respectively.  The  maximum  Q  and  range  occur  at  G^in  “ 

*^min  “  the  value  of  Q  is  a  little  more  than  twice  the  minimum,  and 

the  range  is  indefinitely  large.  These  data  show  that  increased  range 
results  in  an  increase  in  total  heat  input  and  that  the  maximum  total 
heat  input  is  finite  and  occurs  for  an  indefinitely  large  range. 

Figure  15(b)  illustrates  a  family  of  flight  paths,  all  of  which, 
according  to  the  present_analysi5,  have  indefinitely  large  range.  The 
deceleration  begins  at  V  =  2.0  along  a  constant  G  path.  At  a  speci¬ 
fied  velocity  designated  the  flight  path  is  changed  to  a  constant  a 

path,  which  yields  an  indefinitely  large  range.  With  considered  an 

independent  variable,  for  =  2.0  this  case  reduces  to  the  constant  a 

case  previously  considered.  Figure  17  presents  the  total  heat  input 
plotted  against  V^.  Decreasing  reduces  the  total  heat  input.  At 

=  1.3  the  total  heat  input  is  32  percent  more  than  the  value  for  a 
constant  G  of  10  g's,  compared  with  120  percent  more  for  =  2.0, 
which  is  the  constant  a  case. 

To  obtain  control  of  the  range  the  flight  paths  of  figure  15(b), 
like  those  of  figure  15(a),  can  also  be  modified  by  a  constant  G  segment 
through  V  =  1  as  shown  in  figure  15(c).  From  the  discussion  of  the 
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flight  paths  of  figure  15(a),  it  can  be  concluded  that  the  heat  inputs 
for  the  flight  paths  of  figure  15(c)  will  fall  between  the  two  curves  of 
figure  17 . 

Several  coniments  about  the  limitations  of  the  present  range  control 
calculations  may  be  made.  The  composite  flight  paths  considered  in  this 
section  have  mathematical  discontinuities  in  flight-path  angle  at  the 
Junctions  of  the  various  types  of  flight  paths,  so  that  some  transition 
maneuver  is  required.  The  details  of  the  transitions,  however,  are  not 
expected  to  have  a  significant  effect  on  the  total  heat  input  to  the  “i 

vehicle.  c 

c 

,  |- 
The  details  of  the  flight  path  through  V  =  1.0  for  the  constant 

a  flight  plan  and  the  determination  of  a  minimum  permissible  value  of 
both  require  a  more  detailed  study  than  can  be  made  with  the  present 

assumptions . 


Generalizing  the  Results  with  Respect  to  Vehicle 
Gross  Weight  and  Vehicle  Density 

The  numerical  results  for  the  heat  input  have  been  calculated  for 
a  specific  vehicle  with  a  gross  weight  of  10,000  pounds  and  a  vehicle 
density  of  20  pounds  per  cubic  foot.  With  the  assumption  of  a  laminar 
boundary  layer  for  all  vehicles,  these  results  can  be  applied  to  vehicles 
of  different  gross  weight  and  average  density  by  the  following  reasoning. 
Consider  the  results  of  columns  17  and  18  with  the  substitutions  given  at 
the  bottom  of  the  table.  In  the  present  application  the  product 
(Cd/Cr)Cr  is  approximately  For  configurations  like  that  of  fig¬ 
ure  3,  the  drag  coefficient  at  a  =  0,  depends  only  on  the  vehicle 

proportions,  that  is,  on  the  thickness  ratio  d/c  and  sweep  angle  A, 
if  the  contribution  of  friction  is  neglected.  The  contribution  of  fric¬ 
tion  is  small  for  the  low  L/D  confi^ration  considered  in  this  report. 
The  ^/W  for  a  configuration  with  (L/D)inax  =  1  "but  with  any  given 
weight  and  fuselage  density  may  hence  be  estimated  from  the  curves  of 
this  report  by  the  relation: 


X 

curves 


^10,000  20  ^1/6 


CONCLUSIONS 

An  approximate  analysis  of  flight  paths  for  the  descent  from  cir¬ 
cular  and  supercircular  speeds  has  yielded  closed-form  solutions  for  the 
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following  flight-path  characteristics:  altitude,  lift  coefficient,  G 
load,  path  angle,  free-stream  and  local  Reynolds  number'^,  heating  rate, 
heat  input  per  square  foot,  total  heat  input  to  the  vehicle,  range,  and 
time  of  flight.  The  heating  analysis  considers  separately  the  heating 
of  the  leading-edge  region  and  the  windward  side  of  the  vehicle.  The 
six  deceleration  flight  paths  considered  are  constant  values  of  altitude, 
angle  of  attack,  Q  load,  heating  rate,  and  free-stream  Reynolds  number, 
and  a  case  of  modulated  roll.  Limitations  on  altitude,  G  load,  lift 
coefficient,  and  Reynolds  number  have  been  considered.  The  following 
conclusions  may  be  drawn: 

1.  Of  the  flight  paths  considered,  all  except  the  constant-heating- 
rate  path  stayed  within  the  following  physical  limitations:  Altitude 
greater  than  0,  net  acceleration  G  equal  to  or  less  than  10  g's,  and 
lift  coefficient  less  than  0.9  for  flight-  to  circular -velocity  ratios 
equal  to  or  less  than  2.0  and  a  vehicle  maximum  lift-drag  ratio  of  1.0. 
Entry  velocity  ratios  up  to  2.8  are  eierodynamically  possible. 

2 .  Changing  the  type  of  flight  path  or  the  G  level  of  the  flight 
path  has  little  effect  on  the  value  of  the  local  Reynolds  number,  al¬ 
though  the  value  of  the  free-stream  Reynolds  number  may  vary  a  large 
amount.  For  the  present  examples  the  value  of  the  local  Reynolds  number 
is  about  0.7X10®,  a  value  for  which  one  can  be  hopeful  of  achieving  a 
laminar  boundary  layer. 

3.  The  flight  path  with  the  minimum  total  heat  input  is  the  one  with 
the  highest  sdlowable  G  load. 

4.  With  a  laminar  boundary  layer  and  a  maximum  deceleration  of  10 
g's,  entry  velocities  of  twice  circular  and  higher  appear  aerodynamic ally 
and  thermodynamically  feasible  for  a  vehicle  with  a  maximum  lift-drag 
ratio  of  1.0.  About  85  percent  of  the  total  heat  input  to  the  vehicle  is 
to  the  leading-edge  region.  The  ablated  weight  is  estimated  to  be  about 
5  percent  of  the  vehicle  gross  weight  to  decelerate  from  twice  circular 
to  zero  velocity  along  a  constant  10  G  path.  A  constant-angle-of -attack 
path  with  an  initial  G  load  of  4.25  g's  has  a  heat  input  4.9  times 
higher  than  that  of  the  constant  10  G  path. 

5.  For  an  initial  velocity  of  twice  circular,  incorporating  leuiding 
point  control  into  a  constant  10  G  path  by  controlling  the  flight  path 
in  the  vicinity  of  circular  velocity  increases  the  heat  input  about  32 
percent . 


6.  For  a  delta  wing  vehicle  for  which  the  useful  volume  is  the 
volume  in  the  wing,  the  heat  input  per  square  foot  varies  directly  with 
the  one-half  power  of  the  average  fuselage  density  and  inversely  with 
the  one-half  power  of  the  G  load  and  resultant  force  coefficient.  For 
vehicles  of  this  type,  with  the  same  sweep  back,  the  total  heat  input 
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varies  inversely  with  the  one-half  power  of  the  G  load.,  the  one-third 
power  of  the  gross  weight,  and  the  one-sixth  power  of  the  vehicle  average 
density  and  drag  coefficient. 


Lewis  Research  Center 

National  Aeronautics  and  Space  Administration 
Cleveland,  Ohio,  August  7,  1961 
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APPENDIX  A 

SYMBOLS 


reference  area  for  aerodynamic  coefficients,  plan  area  of 
vehicle  (see  fig.  3),  sq,  ft 

projected  frontal  area  of  vehicle  (see  fig.  3),  sq  ft 

span  of  vehicle  (see  fig.  3),  ft 

effective  heat -absorbing  capacity  of  ablation  heat  protection 
system,  Btu/lb 

aerodynamic  dra^  coefficient,  - — ^ - 

ipV^A 


aerodynamic  drag  coefficient  at  zero  lift 
friction  coefficient,  dimensionless 

component  of  aerodynamic  lift  coefficient  in  vertical  direc¬ 
tion, 

ipV-^A 


total  eierodynamic  lift  coefficient 

pressure  coefficient,  2(pg  -  p)/pV^ 

specific  heat  at  constant  pressure,  Btu/( slug) (°R) 

aerodynamic  re sultemt -force  coefficient 
2W 

.  A 

% 

23W 

root  chord  of  vehicle,  ft 

leading-edge  axial  force  coefficient, 

/I  2 

leading-edge  eucial  force/—  pV  A^ 


leading-edge  drag  coefficient, 

ll  2 

leading-edge  drag/—  pV 

leading-edge  lift  coefficient, 
leading-edge  pV^A^ 

leading-edge  normal  force  coefficient, 

ll  2 

leading-edge  normeil  force  pV 

general  wing  chord,  ft 

aerodynamic  drag,  lb 

thickness  of  vehicle  (see  fig.  3),  ft 

see  column  19  of  table  I 

linear  distance  eilong  streamline  from  leading-edge  stagna¬ 
tion  point  to  end  of  leading-edge  region  (beginning  of 
windward  side)  (see  fig.  3),  ft 

net  acceleration,  g's 

acceleration  due  to  gravity,  essentially  the  planet  surface 
value,  ft/ sec 2 

see  column  19  of  table  I 

altitude,  ft 

total  enthalpy  at  edge  of  boundary  layer,  Btu/lb 
wall  enthalpy,  Bbu/lb 

heat-transfer  coefficient,  Btu/(sq  ft)(sec)(°R) 
mechanical  equivalent  of  heat,  778  ft-lb/Btu 
see  eq.  (7a) 

see  columns  11,  17,  18,  and  19  of  table  I 
component  of  aerodynamic  lift  in  vertical  direction 
total  aerodynamic  lift,  lb  force 


E-lOOl 


l/d 

(L/D) 


lift-drag  ratio 
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M 

m 

n 


P 

Pr 


Q 

q 

q 

• 

q 

4 

^l,A=0 

R 

% 

Si 


Re 

Re 
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distance  along  leading  edge  of  wing  (see  fig.  3),  ft 
Mach  number,  dimensionless 
entry  vehicle  gross  mass,  slugs 

exponent  defining  variation  of  altitude  density  with  V  for 
general  flight  path  of  table  I 

static  pressure,  lb  force/sq  ft 

Prandtl  number,  dimensionless 

total  heat  input,  Qj  ®tu 

total  heat  input  to  leading  edge,  Btu 

total  heat  input  to  windward  side,  Btu 

heat  input  per  square  foot ,  Btu/ sq  ft 

dimensionless  function  Q  of  ref.  4,  proportional  to  total 
heat  input  per  unit  area  for  lamlneir  flow 

dimensionless  function  q  of  ref.  4,  proportional  to  heating 
rate  for  laminar  flow 

heating  rate,  Btu/(sq  ft) (sec) 

q^  at  7^0,  maximum  heating  rate  on  swept  leading  edge  (see 
fig.  3  for  A),  Btu/(sq  ft) (sec) 

range,  ft 

leading-edge  radius,  ft 

gas  constant,  1716  ft-lb/(slug) (°R) 

free-stream  Reynolds  number,  dimensionless 

Reynolds  number  at  outer  edge  of  boundary  layer,  dimension¬ 
less 


radius  of  planet,  ft 
Stanton  number,  dimensionless 
total  temperature,  °R 
static  temperature,  °R 
flight  velocity,  ft/sec 

component  of  velocity  normal  to  leading  edge  (see  fig.  3) 

component  of  velocity  normal  to  any  surface  (see  fig.  3) 

ratio  of  flight  to  circuleoc*  velocity,  v/  Vgr,  dimensionless 

dimensionless  velocity  at  which  change  from  constant  G  to 
constant  a  type  of  flight  path  occurs 

entry  vehicle  gross  weight.  Earth  lb 

wing  loading.  Earth  Ib/sq  ft 

ablated  weight,  Earth  lb 

distance  from  leading  edge  in  stream  direction,  ft 
coordinate  systems  on  vehicle  (see  fig.  3) 
angle  of  attack  with  respect  to  free  stream 
cross-flow  angle  of  attack  (see  fig.  3) 

exponent  describing  exponential  variation  of  density  with 

altitude  in  equation  =  e"^^,  ft  (p"^  =  2.35x10^  ft 

PSL 

for  Earth) 

ratio  of  specific  heats,  dimensionless 

angle  of  surface  with  respect  to  free-stream  direction 

emlssivity,  dimensionless 

roll  angle  (see  fig.  l)  (zero  for  lift  in  vertical  plane  di¬ 
rected  upward) 


conductivity,  Btu/(ft) (sec) (°R) 
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A  geometric  leading-edge  sweep  angle  (see  fig.  3) 

Ag  effective  sweep  angle  (see  fig.  3) 

A  leading-edge  radius  coordinate  angle  (see  fig.  3) 

H  viscosity,  slugs/(ft)  (sec)  or  (ITd  force)  (sec)/ sq.  ft 

V  exponent  defining  variation  of  laminar  heating  rate  with  wing 

sweepback  angle 

p  atmospheric  density^  slugs/ cu  ft 

average  entry  vehicle  density.  Earth  Ib/cu  ft 
average  entry  vehicle  density,  slugs/cu  ft 
a  Stefan-Boltzmann  constemt,  0.48X10"^^  Btu/(scL  ft) (sec) (°R^) 

T  time,  sec 

cp  flight-path  angle  with  respect  to  local  horizontal 

(see  fig.  1) 

(0  leading-edge  surface  angle  (see  fig.  3) 

Subscripts;^ 

aw  adiabatic  wall 

e  at  edge  of  boundary  layer 

i  beginning  of  deceleration  flight  path 

I  leading  edge 

meix  maximum 

min  minimum 

n  general  subscript  for  summation  in  appendix  B 

SL  sea  level 

s  stagnation  point 

^Flight  quantities  of  M,  p,  Hi,  Re,  t,  t/T,  V,  f,  p,  and  p  with 
no  subscript  or  superscript  are  free-stream  values. 
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w  windward  side  of  vehicle 

w  wall 

V  vertical  component  of  flight  velocity 

1,2  stations  along  deceleration  flight  path 

Superscript:^ 

*  reference  condition  in  "reference  temperature"  or  "reference 

enthalpy"  method  of  calculating  heat  transfer 


^FlighL  quantities  of  M,  p,  d,  Re,  t,  t/T,  V,  y,  q,  and  p  with 
no  subscript  or  superscript  are  free-stream  values. 
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APPENDIX  B 


VEHICLE  AERODYNAMICS 


The  vehicle  aerodynamics  are  estimated  by  assuming  hypersonic  flow 
and  using  the  approximation  to  the  leading  edge  of  the  vehicle  shown  in 
figure  3. 


The  pressure  coefficient  on  windward  surfaces  is  in  general 


(y  +  l)sin^  =  (y  +  l) 


and  for  a  surface  on  the  leading  edge 


Cp  = 


(r  + 


cos(cu 


a^,)  |cos(a)  -  ttg) 


(Bl) 


(B2) 


In  equation  (B2)  if  Cp  is  positive,  the  surface  is  windward,  if 

negative,  the  surface  is  leeward.  The  pressure  coefficient  on  leeward 
sui’faces  is  assumed  to  be  zero  (rather  than  the  value  given  by  eq.  (B2)). 
From  figure  3 


2  ?  2  2  2 
sin  a  +  cos  a  cos  A  =  1  -  cos  a  sin  A 


(B3) 


and 


a 


c 


eu:c  tan 


tan  g 
cos  A 


(B4) 


Equation  (B2)  hence  becomes 


Cp  =  (y+  l) (l-  cos^a 


sin^A)  COSTCO -arc  tan  ‘g"  cos^cu 


•  arc  tail 


tan  a\ 
cos  A/ 


(B5) 


For  a  leading  edge  composed  of  multiple  flat  surfaces  of  equal  width,  the 
axial  force  coefficient  of  the  leading  edge  using  the  projected  frontal 
area  of  the  vehicle  as  a  reference  area  is 
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'=P,n  % 

1 


COS  (% 


(B6) 


Similarly  the  normal  force  coefficient  for  the  nose  is 

n 

sin  oSjj 


2 


cos  cos  OJijj 

1 


(B7) 


The  contribution  of  the  leading  edge  to  the  drag  and  lift,  from  the  re¬ 
sults  of  equations  (B6)  and  (B7),  is 


Cjj  =  c^  cos  a  +  Cjj  sin  a 


cos  a  -  c^  sin  a 


(B8) 

(B9) 


The  total  lift  and  drag  coefficients  of  the  vehicle  based  on  the  refer¬ 
ence  area  of  the  vehicle  (see  fig.  3)  are  then 


2 

=  (r  +  l)sln  a  cos  o,  +  Cj^  -j- 


Cjj  =  (r  +  l)sin2a  +  Cp  -j- 


(BIO) 

(Bll) 


where  the  first  terms  on  the  right  side  of  the  equations  are  the  contri¬ 
bution  of  the  reference  area  A  to  lift  and  drag  coefficients,  and  the 
second  terms  are  the  contribution  of  the  leading  edge. 

If  one  considers  a  family  of  vehicles  having  the  general  plan  form 

shown  in  figure  3,  that  is,  having  the  same  sweepback  angle  but  various 

values  of  c/d,  or  what  is  equivalent,  various  values  of  (L/D)  ,  then 

m&x 

the  following  relation  may  be  derived  for  low  emgles  of  attack: 

Af 

~  ^D,o  =  =D,o  X 


E-lOOl 


E-lOOl 


35 


where  o  constant  for  a  given  sveephaok  angle. 

(El) 


A  « 


and 


Af.  =  bd 


From  equation 


so  that 


^  I  =  2°D,o  =  Constant 


(B12) 


APPENDIX  C 


WINDWARD-SIDE  HEATING  ANALYSIS 


A  simple  analytic  expression  for  the  laminar  heat  transfer  to  the 
windward  side  of  the  vehicle  at  high  velocities  is  developed  herein. 

Off  the  stagnation  point  for  moderate  angles  of  attack  the  heat 
transfer  to  the  v/indward  side  is  analyzed  using  an  approximation  to  the 
"reference  temperature"  (ref.  9).  The  relations  required  are: 

4  = 


St  = 


^  e  p 


St  =  -y  (Pr*)*^/^ 


(C2) 

(C3) 


Pr* 


(C4) 

(C5) 

(C6) 


t*  =  tg  +  0.5(t^  -  tg)  +  0.22(ta„  -  tg) 


(C7) 


Assuming  t^  «  t^^  and  combining  equations  (Cl)  to  (C5)  yield 


<1  = 


0.166  2 
jpy*0.167l  X 


(C8) 


The  three  terms  in  the  numerator  of  the  square  root  term  of  equation 
(C8)  may  be  evaluated  as  follows,  as  an  approximation  of  Southerland's 
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viscosity  equation  for  high  temperatures  (ref.  17,  eq.  (A3)) 


=  2.27x10"®  2.27X10"®  t*V2 


(C9) 


The  reference  density  is  given  by 


*  t 

p  =  p  T  ^ 


(CIO) 


At  moderate  angles  of  attack,  by  Newtoniem  aerodynamics. 


Vg  =  V  cos  a 


(Cll) 


The  terms  t*  and  Pg/p  in  equations  (C9)  and  (CIO)  must  now  be 
evaluated.  In  equation  (C7)  for  low  and  moderate  angles  of  attack,  tg 
and  (t^  -  tg)  are  small  ccmpeired  with  (t^^  -  tg)  and  t^^,  so  that  equa¬ 
tion  (C7),  combined  with  equation  (C5),  beccmes 


f  .  0.22  t  =  0-UVf?^V 

P 


(C7a) 


From  hypersonic  aerodynamics 


Cp  e  =  {—  -  1)  — ^  =  (r  +  l)sin^a 

VP  /  Im® 


(C13) 


where  at  supercircular  speed  (M  >  25)  and  for  moderate  angles  of  attack 
Pe/P  I'®'  angles  of  attack  where  the  assumption  produces 

maximum  error,  the  heat  input  to  the  windward  side  of  the  vehicle  is  gen¬ 
erally  small  compared  with  the  heat  input  to  the  leading  edge  (see  fig. 
14) .  The  error  in  estimating  the  total  heat  input  to  the  vehicle  can 
thus  still  be  small  at  small  angles  of  attack.  Then,  since 


equation  (C13)  becomes 


m2  V 
Yi^t 


°e  ,  .  r  +  1  V  , 

—  =  1  +  —  ST  sin'^a 

p  2  ftt 


r  +  .P.  ri  Sin2a 

2  ftt 


(C14) 
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Making  the  appropriate  substitution  In  equation  (C8)  gives 


0.166 


fe. 27X10“® +  1  pV^  cos  a  sin^a  7t2 


j0.75pr*0.292'y 


gr  (CIS) 


The  following  values  are  assumed  to  be  constants : 

Pr*  »  0.71 

r  "  1.2 

ft  =  1716  ft-lb/(slug)(OR) 

Cp  =■  j'Cy^  i)  “13.2  Btu/(slug)(°R) 

J  =  778  ft-lb/Btu 

From  the  hypersonic  aerodynamics  of  the  configuration,  equation 
(Bll),  and  equation  (6) 


A  cos  6  A 


^a  r  +  1  „2  .  2 

-T-  »  ■* — s-  pV  sin  a  cos  a 


(C16) 


where 


K  =  |1  +  — 
"  \  A(r 


Vl  Y 

+  l)sln^a  cos  a/ 


so  that  equation  (ClS)  becomes,  with  the  assigned  constants. 


(C17) 


(C18) 


The  absolute  value  signs  are  used  to  Indicate  that  the  heating  Is  real 
and  positive  irrespective  of  the  sign  of  L  and  cos  0.  By  use  of  the 

vedue  of  -  ^  given  In  column  6  of  table  I,  one  may  solve  for  L  In 
g  dr  ^  -w 

equation  (4b)  and  substitute  In  equation  (C18)  to  arrive  at  the  results 
given  in  column  10  of  table  I.  For  example,  for  the  constemt  Re  case 
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with  0=0 

^  =  1.50xl0-\g3/2r^^l  _  v2  V  v2»10.2  ^  -  %e  ^ 

(C19) 

2 

The  term  k^,  is  the  ratio  of  the  lift  produced  hy  the  vehicle  wind¬ 
ward  side  to  the  sum  of  the  windward-side  eind  leading-edge  lift.  Approx¬ 
imately,  for  angles  of  attack  from  0°  to  15° 

f  Af  dCj^  2^  Vl/2 

A  d(8in  a)  (y  +  l)sln  aj 

where  dcL/d(sin  a)  is  approximately  constant  and  depends  on  the  sweep 

angle.  For  the  present  configuration  with  A  =  60°  and 
dcjyd(sin  a)  =  0.88,  A^/A  =  0.569,  a  «  5.3°,  and  k^  =  0.535. 

The  term  k^,  in  genered,  varies  over  the  flight  path.  For  super¬ 
circular  entry  velocities  a  conservative  estimate  of  the  heat  input  is 
obtained  by  evaluating  k^  at  the  maximum  velocity.  This  is  done  in 

the  present  calculations.  The  value  of  kc  approaches  unity  and  con¬ 
stancy  for  small  values  of  A^/A  or  high  values  of  (L/d)^^^ . 


40 


APPENDIX  D 


METHOD  USED  FOR  CALCULATING  RATIO  OF  LOCAL 
TO  PREE-STREAM  REYNOLDS  NUMBER 

The  procedure  for  calculating  the  local  Reynolds  number  at  super¬ 
sonic  speeds  for  cones  and  wedges  with  unswept  edges  is  presented  in 
reference  18.  A  relation  for  calculating  the  local  Reynolds  number  at 
hypersonic  speeds  for  the  case  of  a  flat  plate  at  angle  of  attack  with  a 
swept  blunt  leading  edge  is  developed  herein.  The  nomenclature  is  illus¬ 
trated  in  figures  3  and  18. 

At  a  station  sufficiently  far  downstream  from  the  leading, edge, 
where  the  surface  static  pressure  for  a  blunt  body  is  approximately  that 
for  a  nonblunt  body,  the  ratio  of  local  to  free-stream  Reynolds  number 
per  foot  or  for  the  same  length  may  be  written 


Re^  _  Pe  u 
W  “  p  V  |ig 


(Dl) 


In  equation  (Dl)  the  velocity  ratio  is  taken  as  the  Newtonian  value. 


cos  a 


(D2) 


(C9) 


The  viscosity  approximation  is  the  same  as  that  used  in  equation 

> 


e 


2.27X10"® 


(D3) 


where  for  this  case 


t 


(D4) 


where 


V2  cos^A  cos^A 

t  _  f  u.  _ L  %  —  — .  ^ 

^  ^  2JCp  2JCp 


(D5) 
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Hence  equation  (D3)  'becanes 


H-  =  2.27X10' 


The  term  Pg/p  in  equation  (Dl)  may  he  written 


^  _  Ps  Pe 
P  "  P  Ps 


(D6) 


(D7) 


and  equation  (Dl)  then  heccoies 
Re 


'I 


M  cos  aV2JC„  — 

F  £_ 


■(^) 

c '  s' 


¥ 


2.27XL0-8  VVgr  cos  A 


The  ratio  Pg/Ps  written  by  use  of  the  perfect  gas  law  as 


Pe  /Pe\ 

Ps  '  W  '  \P/ 


Hypersonic  aerodynamics  gives  (see  also  eq.  (C14)) 


(D8) 


(D9) 


Pe  T  j.  r  +  1  ,  2 

T  °  ^  2~  w  “ 


(DIO) 


and 


^  =  1  +  y  ^  cosSa  *  IT  +  1  YI  ggs2A 
p  2  Et  c  2  fit  c 


(Dll) 


where  the  approximation  is  valid  for  sweep  angles  less  than  about  87°. 
From  figure  3  (see  also  eq.  (B3)) 


cos^Ag  =  (l  -  sin^A  cos^a) 


(D12) 


The  density  ratio  across  the  shock  at  the  cross-flow  stagnation  point  is 


^  =  r  +  ^ 

P  r  -  1 


(D13) 
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With  the  appropriate  substitution,  equation  (D8)  gives  finally 


Re, 


[i  cos  aV2JCp^|^^ 


2. 27X10-®  VV^(1  -  sin^A  oos^a)^/^ 


2ftt 


,2r 


(r  +  l)^gr 


+  sin®a 


1  -  sin^A  cos^a 


J  (D14) 


The  following  vsiLues  were  assumed  in  the  present  calculations: 

Y  =  1.2 

ft  =  1716  ft-lb/(slug)(^) 

Cp  =  2.)  “  Btu/(slug)(°R) 

t  =  500°  R 

pi  =  0.37X10"®  (lb  force)  (sec)/ sq  ft 
J  =  778  ft-lb/Btu 

The  term  V  is  the  independent  variable,  and  a  is  determined  by  the 
lift  coefficient  (column  7  of  table  l)  and  the  vehicle  aerodynamics. 

This  celculatlon  takes  no  detailed  account  of  the  boundary-layer 
history  on  the  leading  edge  of  the  vehicle. 

Typical  curves  of  the  ratio  of  local  to  free-stream  Reynolds  number 
are  shown  in  figure  18  as  a  function  of  angle  of  attack  for  several 
values  of  V.  Both  V  and  angle  of  attack  have  a  strong  effect  on  the 
Reynolds  number  ratio. 
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APPENDIX  E 


HfTECaiATION  OF  HEAT  INPUT  OVER  VEHICLE 

The  vehicle  geometry  over  which  the  integration  of  the  local  heat 
per  square  foot  q  is  carried  out  is  shown  in  figure  3  together  with 
the  symbols  and  coordinate  systems  used.  The  constant  Re  case  is  used 
as  an  example. 


Vehicle  Geometry 

The  following  approximate  relations  can  be  derived  from  the  vehicle 
geometry : 
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Integration  Over  Leading-Edge  Region 


The  integration  over  the  vehicle  is  divided  into  two  parts,  that 
over  leading-edge  regions,  and  that  over  the  windward  side  of  the  ve¬ 
hicle.  To  illustrate  the  integration  over  the  leading-edge  regions, 

may  be  taken  from  column  15  of  table  I  (hence  the  subscript  15  on  k)  or 
the  integration  of  equation  (37)  and  written  as 


,v/2 


q^  =  kj^gfl  -  sin^A  cos^a]  cos  A 
where  for  constant  Re 

r^.5  _ 

H  1  L  3.5 


(E2) 


ki5  =  12.1X10^ 


Vkv) 


The  general  form  of  the  integration  desired  is  (see  fig.  3) 

n  /•  n/2 


•'o  Mrt 


(jt/2-ac) 


(E3) 


Carrying  out  the  interior  integration  and  using  equations  (El)  and  (E2) 
yield 


/ 


'/2 


q,i  =  k^gd  /  (l  -  sin^A  cos^aj  (l  +  cos  a,g)dl 


(E4) 


where  from  figiire  3 


tan  a 

a„  =  arc  tan  - r 

c  cos  A 


Equation  (E4)  may  be  written 


/- 


cos^A  dZ 


(E5) 


(E6) 
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where 


|l  -  sin^A  cos^a|  ^ 


1  +  cos  arc  tan 


(tan  a\) 

COB  A/f 


2  cos^A 


(E7) 


For  a  =  0,  kj  =  1.0;  for  a  =  12°,  v  =  1.5,  and  A  =  60°,  k^  =  1.05; 

thus,  for  the  integrations  over  the  flight  paths  for  which  the  angle  of 
attack  is  small,  kj  is  assumed  to  be  unity.  For  the  constant  a  case, 

which  is  used  at  higher  a,  kj  is  a  constant  and  may  be  evaluated.  With 
kj  a  constant  and  by  use  of  equation  (El),  equation  (E6)  can  be  written 


(v-1)a 


/b/2 

db 


so  that  finally  integration  gives 

Qj  =  ^  kiskjbd  cos('^“^^A  =  ^  k25k2A|.  cos^^'^^A 


(E8) 


(E9) 


Writing  this  result  in  terms  of  the  gross  weight  and  wing  loading,  by 
using  equation  (El) ,  yields 


^  _  kigk^  cob(v-1)a  /  w/A 
W  ~  \|W  tan  A/ 


(ElO) 


or 

^  _  k^gk^  cos(^-1)a  / 
m  ”  \pi  tan  A/ 

and  substituting  for  k^^g  (eq.  (E2))  and  for  W/Rj^A  from  equation  (El) 
results  in 
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or 


Jz(.^.3JA _ .\ 

-1)^  \PjHi  tan  hj 


1/2 


cos^'^ 


Comparing  equations  (E2)  emd  (Ell)  yields 


[^2.5  yS.sl 


W 


<1 


I 


( 1  -  sin^A  cos^a) 


v/2 


cos  X 


•V2kj  cos^'^'^^A 
/pjW  tan  A\^/^ 
V  w/A  ; 


(E12) 


Equation  (E12)  results  from  an  integration  over  the  vehicle  end  hence 
is  true  for  all  flight  paths.  It  is  listed  in  column  17  of  table  I  and 
a  similar  relation  is  given  in  table  II. 


Integration  Over  Windward  Side  of  Airplane 

For  this  case  q^,  may  be  taken  from  column  16  of  table  I  (hence 

the  subscript  16  on  k)  or  the  Integration  of  equation  (30)  and  may  be 
written  as 


^  * 


(E13) 


where,  for  example,  for  the  constant  Re  case  for  V  >  1, 


The  integration  desired  is  ther.  (see  fig.  3) 


Vs  "  T)  tan  A 


(E14) 
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The  Integration  results  in 

8 


^  “  3  ^16 
Equation  (EIS)  may  be  written 


(f  4-  b 

tan  A  "  tan  A  “  4 


(E15) 


8,  ,  c3/2 


^  =  3  1^16^  tan  A 
where,  using  equation  (El), 


(E16) 


f\3/2  /f\3/2  ^  3/f 


(E17) 


In  terms  of  gross  weight  and  wing  loading  and  by  use  of  relations  from 
equation  (El), 


^  8  .  .  _  1 

w  ■  3  —374 


(E18) 


(W  tan  A)^/^ 


Substituting  for  (eq.  (E13))  gives 


(V^  -  1  +  ^ 


Lv) 

5%  {l  ^ 


(E19) 


Comparing  equations  (E13)  and  (E19)  gives 

‘Iw  8  ^ 


^  (I  W  tan  a) 


174 


(E20) 


Equation  (E20)  holds  for  sll  flight  paths  and  is  given  in  column  18  of 
table  I.  A  similar  result  is  given  in  table  II. 
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Prom  section  AA  of  figxire  3  the  ratio  f/c  may  be  written  approxi¬ 
mately  for  small  angles  of  attack  as 

and  by  use  of  equation  (E5) 


f 

c 


f-p 

4  c  L 


cos  A 


4 

—  arc  tan 
n 


tan  g 
cos  A 


(E22) 


For  the  configuration  used  as  an  example  in  this  report  to  generate 
numerical  values,  %  =  2  .06,  c  =  14.5,  and  A  »  60°.  For  a  »  0°  and 

5.3°,  f/c  is  0.34  and  0.31,  respectively.  The  corresponding  values  of 
(eq.  (E17))  are  0.48  and  0.49,  respectively.  Note  that  ky  does  not 
vary  significantly  for  small  changes  in  a,  and  hence  it  is  assumed  to 
be  constant.  A  value  of  0.49  was  used  for  the  calculations  presented. 

The  most  desirable  form  of  seme  of  the  terms  in  equations  (Ell)  and 
(E19)  depends  on  what  information  is  known.  The  following  useful  alter¬ 
nate  forms  may  be  derived  if  the  parameter  W/A  is  eliminated  by  use  of 
one  of  the  equations  (El): 


(E23) 

(E24)' 


The  term  c/d  depends  on  the  vehicle  sweep  angle  and  maximum  lift -drag 
ratio.  For  the  present  examples  with  A  e  60°,  c/d  =  3.51  and  21.7  for 
maximum  lift-drag  ratios  of  1  and  2,  respectively.  Still  another  inter¬ 
esting  form  of  the  results  of  columns  17  and  18  mey  be  written  using  the 
results  of  equations  (E23) ,  (E24),  and  (B13) .  For  vehicle  operation  at 
low  engles  of  attack  where  ^  ~  ~  which  is  approximately  the 

same  assumption  made  to  arrive  at  the  results  of  columns  17  and  18, 


-  cJ/V/V/®(2=b  tan 


(E25) 


TOOT- a 


E-lOOl 
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(«  W  tan  -  cj/V/ipys  (E26) 

where  Cjj  is  constant  for  a  given  sweepback  angle. 

An  example  of  the  significance  of  this  result  is  the  following.  For 
a  family  of  vehicles  having  the  same  sweepback  euigle,  gross  weighty  and 
fuselage  density,  the  total  heat  input  to  the  vehicle  leading  edge  varies 
as 


(E27) 


(by  the  results  of  column  17  and  eq.  (E25)  for  kj  =  Constant) .  On  the 

other  hand  (by  column  15)  the  heat  input  per  square  foot  to  the  leading 
edge  varies  as 


(E28) 


Thus  the  change  in  the  total  heat  input  to  the  vehicle  is  much  less 
sensitive  to  changes  in  Cp  (or  (L/D)^^^^)  than  is  the  heat  input  per 
square  foot. 
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APPENDIX  F 


CASE  OF  CONSTANT  ANGLE  OF  ATTACK 


The  determination  of  the  term  —  ^  is  more  difficult  for  the  case 

g  dr 

of  constant  angle  of  attack  than  for  the  other  flight  paths.  However, 
the  solution  of  this  case  permits  a  check  on  the  accuracy  of  the  results 
of  the  present  analysis  hy  comparison  with  the  numerical  integrations  of 
references  4  and  13.  The  equations  to  be  manipulated  are  repeated  here 


1  ^  ^ 

i  dT  =  -  ^  C^ 

0^=  1  .  v2 

^  g  dT 

CRpV^Agr 
°  “  - 2W - 

CjjgVA  ^ 

^  “  2PW 

Using  the  mathematical  identity 

V  ^  _  V  ^  ^ 
g  dr  8  dv 

and  equation  (3)  permits  equation  (4a)  to  be  written 


(3) 

(4a) 


(5) 

(19) 


(FI) 


Solving  equation  (5)  for  p,  by  use  of  (F2)  for  G,  yields 


(F2) 


-  ‘=p(p  -  i)  f 


(F3) 


E-lOOl 
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where 


C  H— 2W_ 
P  Cl 

CRAgr  ^ 

Differentiating  equation  (19)  gives 


where 

c 

<P  “  23W 


(F4) 


(F5) 


(P6) 


The  term  dq)/dV  can  be  eliminated  between  equations  (F3)  and  (F5)  to 
yield  a  single  differential  equation  fca*  which  a  solution  must  be  found. 
Herein,  however,  a  solution  to  equations  (F3)  and  (F5)  was  found  by 
successive  approximations.  For  a  first  approximation  the  flight-path 


angle  is  assumed  to  be  a  constant 


^  »  oV  Then  from  equation 
dV  / 


(P3) 


and 


dV^ 


and  equation  (F5)  becomes 


dtp  ^  _  4 

dV  “  prV^  ^ 


(F7) 


(F8) 


(F9) 


The  value  for  dqp/dV  given  by  equation  (F9)  may  now  be  introduced  into 
equation  (F3)  to  obtain  a  second  approximation.  Thus, 


(FIO) 
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where 


Pr(L/D)2 


(Fll) 


Repeating  the  steps  of  equations  (F7)  to  (F9)  results  in 


Higher  order  approximations  may  be  obtained.  If  this  degree  of  approxi¬ 
mation  is  accepted,  equation  (F2)  for  G  becomes 


G  = 


(F15) 


where 


F 


1  + 


k 

72  +  k 


1 


+ 


(1  -  7^)k/  2k 

72  +  k  \72  +  k 


■1 


(F16) 


With  G  known,  p  can  be  determined  from  equation  (5),  and  this  is  the 
result  given  in  table  I. 


From  equations  (19)  and  (F12) 


where 


2H 


1^(1  - 
72  +  k  J 


(F17) 


(F18) 


E-lOOl 


Prom  equation  (Pl),  'by  use  of  equations  (3)  and  (P14), 


which  Is  the  result  given  in  column  6  of  table  I. 


APPENDIX  G 


COMPARISON  WITH  OTHER  ANALYSES 


The  results  of  reference  4  are  essentially  for  a  three-dimensional 
stagnation  point,  or  J  =  1  In  equation  (7a).  The  result  of  column  10 
of  table  I  for  constant  angle  of  attack  for  a  three-dimensional  stagna¬ 
tion  point  can  be  written  in  the  modified  form 


4  -  V2  14.9->/i 


V^N^  Vl 


-/ITd  Vf 


(01) 


Prom  equation  (36)  of  reference  4  for  Earth,  and  in  the  nomenclature  of 
this  report. 


Substituting  equation  (Gl)  in  (G2)  gives 


(02) 


4  -  0.202 


(03) 


The  equation  from  column  15  of  table  I  for  a  three-dimensional  nose  is, 
in  modified  form. 


1  -  V2  12. 05X10=  Vi 


■Tv, 


arc  sin  V  -  vVl  - 


(04) 


Prom  equation  (39a)  of  reference  4,  and  in  the  nomenclature  of  this 
report. 


q 


3. 


15,900 


(05) 
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Evaluating  (04)  for  »  0.99  and  V2  "  0, 
4,  and  substituting  in  equation  (05)  yield 


the  values  used  in  reference 


«1 


3.91 


(G6) 


As  a  point  of  interest,  the  total  heat  input  per  square  foot  q  or  q 
is  quite  sensitive  to  the  initi^  velocity  for  this  case.  For  example, 
evaluating  equation  (04)  from  =  1.0  to  V2  ■  0  gives  for  equation 

(G6) 


q=4.76^^  (07) 

which  is  an  Increase  of  22  percent  over  the  case  for  =  0.99.  In  the 

£resent  calculation  the  term  F  is  evaluated  at  the  highest  value  of 
V.  If  the  flight  path_i8  divided  into  several  parts,  F  is  evaluated  at 
the  highest  value  of  V  in  each  part. 

Reference  4  also  presents  the  horizontal  component  of  the  0  load, 
0|j.  For  small  flight-path  angles,  an  assumption  already  made  in  the 

present  analysis. 


Oh  «  0  ^  (08) 

A 

V.'ith  the  result  given  for  0  in  column  4  of  table  I  for  the  constant  a 
case. 


(09) 


Part  of  the  difference  in  heating  rate  between  the  present  analysis 
and  reference  4  is  due  to  the  use  of  a  slightly  different  value  of  the 
constant  in  equation  (7a) . 
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0  20  40  60  80  100 

Angle  of  attack,  a,  deg 

(a)  Lift  coefficient. 


Figure  2.  -  Aerodynamic  characteristics. 


20  40  60  80  100 

Angle  of  attack,  a,  deg 


(b)  Drag  coefficient. 

Figure  2.  -  Continued.  Aerodynamic  characteristics. 
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0  20  40  60  80  100 

Angle  of  attack^  a,  deg 

(c)  Resultant -force  coefficient. 

Figure  2.  -  Continued.  Aerodynamic  characteristics. 
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(d)  Ratio  of  lift  to  resultant  force. 


Figure  2.  -  Continued.  Aerodynamic  characteristics. 
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Figure  2. 


40  60 

Angle  of  attack,  a,  deg 

(f)  Lift -drag  ratio. 

Concluded.  Aerodynamic  characteristics. 


Figure  3.  -  Vehicle  geometry,  coordinate  systems,  and  flow  angles. 
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Reference  area  for  aerodynamic 
coefficients  and  wing  loading; 
also  used  for  volume  estimation 


Figure  3, 


-  Concluded.  Vehicle  geometry,  coordinate  systems,  and  flow  angles. 


Hondlmensional  total  heat  Hondlmensional  heating  rate  at  Horizontal  component  of 

input  at  spherical  spherical  stagnation  point,  G  load,  Og,  g' s 

stagnation  point,  4 
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Present  analysis 
Ref.  15,  fig.  4 
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Fllf^ht-  to  circular-velocity  ratio,  V 


Present  analysis 
Ref.  13,  fig.  7 


— Present  analysis 
-  Ref.  13,  fig.  17 


Vehicle  angle  of  attack,  a,  deg 

Figure  1’.  -  Comparison  of  present  analysis  for  windward-side  heating 
with  results  of  reference  13,  Initial  velocity  ratio,  l.Oj  wall 
static  temperature,  0®  R>  constant -angle-of -attack  flight  path. 
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Time,  T,  min 


Figure  6.  -  Conetant-angle-of-attack  flight  path  for  Earth  atmoa- 
phere.  Angle  of  attack,  5.3°(  Initial  net  acceleration,  10  g'a; 
Initial  velocity  ratio,  2.0;  maximum  lift-drag  ratio,  1.0;  entry 
vehicle  gross  weight,  10,000  pounds;  average  entry  vehicle  den¬ 
sity,  20  pounds  per  cubic  foot;  geometric  leading-edge  sweep 
angle,  60°. 
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Ftatlo  of  flight  to  circular  velocity,  V 


Figure  6.  -  Concluded.  Constant-angle-of -attack  flight  path  for  Barth 
atmosphere.  Angle  of  attack,  5.3°j  Initial  net  acceleration,  10  g'sj 
initial  velocity  ratio,  2.0|  maximum  lift-drag  ratio,  l.Oj  entry  vehicle 
gross  weight,  10,000  pounds;  average  entry  vehicle  density,  20  pounds  per 
cubic  foot;  geometric  leading-edge  sweep  angle,  60°. 


Ratio  of  flight  to  clro  ilar  vploclty,  V 

O.IO  2.70  .'lO  0  -.54  -.90 

Time,  T,  min 


Fig  ii'e  I.  -  Ci'ii3taiit-aitlt  jde  fllgta  jjath  frr  Eai’th  atmosphere. 
Altltide,  102,000  feet;.  Initial  net  .acceleration,  10  g'sj 
lrrltl.al  velr.clly  ratio,  2.0j  maxl'i.ui’.  I1ft-drng  ratio,  1.0; 
entry  vehicle  gross  weight,  10,000  ptuiids;  average  entry  vehicle 
density,  20  pfinnds  per  cthlc  fr.r.t;  ,^ef.t»elrlo  leadlnr-edge  sweep 
angle,  60°. 


Time,  T,  min 

Plpure  0.  -  Constant  free-atream  Reynolds  number  flight  i»th  for  Earth  atmos 
phero.  Free-atream  Reynolds  number,  7.3X10®J  Initial  net  acceleration, 

10  p's)  initial  velocity  ratio,  2.0;  maximum  lift-drag  ratio,  1.0;  entry 
vehicle  grosa  weight,  10,000  pounds;  average  entry  vehicle  density,  20 
pounds  per  cubic  foot;  geometric  leading-edge  sweep  angle,  60°. 
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Ratio  of  flight  to  circular  velocity,  V 


Figure  8.  -  Concluded.  Constant  free-stream  Reynolds  numher  flight  path  for 
Earth  atmosphere.  Free-stream  Reynolds  number.  7.3X10  j  Initial  net  accel¬ 
eration,  10  g'sj  initial  velocity  ratio,  n.O}  maximum  lift-drag  ratio.  1.0; 
entry  vehicle  p.ross  weight,  10,000  poiuids;  averap.e  entry  vehicle  density, 

20  pounds  per  cubic  foot;  geoirietric  leading-edge  sweep  angle, 
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Ratio  of  fli/;^ht  to  circular  velocity,  V 
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Time,  T,  min 


Figure  'J.  -  Constant-net-acceXeratlon  path  for  Eartli  atmosphere.  Net 

acceleration,  10  r/s;  Initial  velocity  ratio,  2.0i  maximiam  llft-draf’;  ratio, 
1,0;  entry  vehicle  f^ross  weight,  10,000  pounds;  averui'^e  entry  vehicle  den¬ 
sity.  20  pounds  per  cubic  foot;  n^ometrlc  leadinf^edpe  sweep  anp:le.  60°. 
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Ratio  of  night  to  circular  velocity,  V 

*  Figure  9.  -  Concluded.  Constant-net-acceleratlon  flight  path  for 

Earth  atmosphere.  Net  acceleration,  10  g's;  Initial  velocity 
ratio,  2.0j  maximum  lift-drag  ratio,  1.0;  entry  vehicle  gross 
weight,  10,000  pounds;  average  entry  vehicle  density,  20  pounds 
per  cubic  foot;  geometric  leading-edge  sweep  angle,  60°. 
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Ratio  of  night  to  circular  velocity,  V 

Figure  10.  -  Constant-heatlng-rate  flight  path  for  Earth  atmosphere.  Maxi 
mum  leading-edge  heating  rate,  900  Btu  per  square  foot  per  second;  Ini¬ 
tial  net  acceleration,  10  g's;  Initial  velocity  ratio,  2.0;  maximum  lift 
drag  ratio,  1.0;  entry  vehicle  gross  weight,  10,000  pounds;  average 
entry  vehicle  density,  20  pounds  per  cubic  foot;  geometric  leading-edge 
sweep  angle,  60°. 
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Figure  10.  -  Concluded.  Constant-heatlng-rate  flight  path  for 
Earth  atmosphere.  Maximum  leading-edge  heating  rate,  900  Btu 
per  square  foot  per  second;  Initial  net  acceleration,  10  g's; 
Initial  velocity  ratio,  2.0;  maximum  lift-drag  ratio,  1.0; 
entry  vehicle  gross  weight,  10,000  pounds;  average  entry 
vehicle  density,  20  pounds  per  cubic  foot;  geometric  leading- 
edge  sweep  angle,  60°. 


Figure  11.  -  Modulated  roll  flight  path  for  Earth  atmosphere.  Initial  roll 
angle,  0;  Initial  net  acceleration,  10  g's;  initial  velocity  ratio.  2.0j 
maximum  llft-d-^ag  ratio,  l.Oj  entry  vehicle  gross  weig'ut,  10.000  pounds; 
average  entry  vehicle  density,  20  pounds  per  cubic  foot;  geometric  leading 
edge  sweep  angle,  60°. 
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Figure  12.  -  Total  heat  input  at  Earth  for  flight  paths  with 
net  acceleration  of  10  g's  at  velocity  ratio  of  2.0.  Maxi¬ 
mum  lift-drag  ratio,  1.0)  entry  vehicle  gross  weight,  10,000 
pounds;  average  entry  vehicle  density,  20  pounds  per  cubic 
foot)  geometric  leading-edge  sweep  angle,  60°. 
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Figure  13.  -  Total  heat  input  at  Earth  for  flight  path  with  Initial 
net  acceleration  of  10  g's.  Maximum  lift-drag  ratio,  1.0;  entry 
vehicle  gross  weight,  10,000  pounds;  average  entry  vehicle  den¬ 
sity,  20  pounds  per  cubic  foot;  geometric  leading-edge  sweep 
angle,  60^\ 
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Figure  16.  -  Total  heat  input  at  Earth  for  modified  constant -angle- 
of -attack  flight  path  (fig.  15(a)).  Initial  net  acceleration, 

10  g*s;  initial  velocity  ratio,  2.0;  maximum  llft-drsig  ratio,  1.0; 
entry  vehicle  gross  weight,  10,000  pounds;  average  entry  vehicle 
density,  20;  geometric  leading-edge  sweep  angle,  60°. 
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Figure  17.  -  Total  heat  input  at  Earth  for  modified 
constant-net -acceleration  flight  path  (fig.  15(b)) 
and  infinite  range.  Initial  net  acceleration,  10  g’sj 
initial  velocity  ratio,  2.0;  maximum  lift-drag  ratio, 
1.0;  entry  vehicle  gross  weight,  10,000  pounds;  aver¬ 
age  entry  vehicle  density,  20;  geometric  leading-edge 
sweep  angle,  60°. 
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Flgvire  18.  -  Ratio  of  local  to  free-stream  Reynolds  number  for 
swept  blunt  leading  edge  for  Earth  atmosphere.  Leading-edge 
sweep  angle,  60)  specific-heat  ratio,  1.2. 
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